arXiv:hep-ph/0503272v4 27 Sep 2009 


Pion Superfluidity and Meson Properties at Finite Isospin Density 


Lianyi He, Meng Jin, and Pengfei Zhuang 
Physics Department, Tsinghua University, Beijing 100084, China 


We investigate pion superfluidity and its effect on meson 
properties and equation of state at finite temperature and 
isospin and baryon densities in the frame of standard flavor 
SU(2) NJL model. In mean field approximation to quarks and 
random phase approximation to mesons, the critical isospin 
chemical potential for pion superfluidity is exactly the pion 
mass in the vacuum, and corresponding to the isospin sym¬ 
metry spontaneous breaking, there is in the pion superfluid¬ 
ity phase a Goldstone mode which is the linear combination 
of the normal sigma and charged pion modes. We calculate 
numerically the gap equations for the chiral and pion con¬ 
densates, the phase diagrams, the meson spectra, and the 
equation of state, and compare them with that obtained in 
other effective models. The competitions between pion super¬ 
fluidity and color superconductivity at finite baryon density 
and between pion and kaon superfluidity at finite strangeness 
density in flavor SU(3) NJL model are briefly discussed. 
PACS : ll.lO.Wx, 12.38. - t, 25.75.Nq 

I. INTRODUCTION 

It is generally believed that there exists a rich phase 
structure of Quantum Chromodynamics (QCD) at finite 
temperature and baryon density, for instance, the decon¬ 
finement process from hadron gas to quark-gluon plasma, 
the transition from chiral symmetry breaking to the sym¬ 
metry restoration [1], and the color superconductivity [2] 
at low temperature but high baryon density. Recently, 
the study on the QCD phase structure is extended to 
finite isospin density. The physical motivation to study 
QCD at finite isospin density and the corresponding pion 
superfluidity is related to the investigation of compact 
stars, isospin asymmetric nuclear matter and heavy ion 
collisions at intermediate energies. In early studies on 
dense nuclear matter and compact stars, it has been sug¬ 
gested that charged pions and even kaons are condensed 
at sufficiently high density [3-6]. 

While the perturbation theory of QCD can well de¬ 
scribe the properties of the new phases at high tem¬ 
perature and/or high density, the study on the phase 
structure at moderate (baryon or/and isospin) density 
depends on lattice QCD calculation and effective models 
with QCD symmetries. While there is not yet precise 
lattice result at finite baryon density due to the Fermion 
sign problem [7], it is in principle no problem to do lat¬ 
tice simulation at finite isospin density [8]. Recently it is 
found that [9-11] there is a phase transition from normal 
phase to pion superfluidity at a critical isospin chemical 
potential which is about the pion mass in the vacuum, 


/Tj ~ ttItt. The QCD phase structure at finite isospin 
density is also investigated in many low energy effec¬ 
tive models, such as chiral perturbation theory [8,12-15], 
Nambu-Jona-Lasinio (NJL) model [16-20], random ma¬ 
trix method [21,22], ladder QCD [23], and strong cou¬ 
pling lattice QCD [24]. 

One of the models that enables us to see directly how 
the dynamic mechanisms of chiral symmetry breaking 
and restoration operate is the NJL model [25] applied to 
quarks [26-29]. The chiral phase transition line [26-32] in 
the temperature and baryon chemical potential {T — ^b) 
plane calculated in the model is very close to the one 
obtained with lattice QCD. Recently, this model is also 
used to investigate the color superconductivity at moder¬ 
ate baryon density [33-40]. It is nature to extend the NJL 
model to finite isospin chemical potential [16,17]. In the 
frame of general NJL Lagrangian with 1 / 4 ( 1 ) breaking 
term [19], it is analytically proved [20] at quark level that 
the critical isospin chemical potential for pion superfluid¬ 
ity is exactly the pion mass in the vacuum, /rj = m^n-, and 
this relation is independent of the model parameters, reg¬ 
ularization scheme, 1/^(1) breaking coupling, and color 
degrees of freedom N^. 

Unlike the effective models at hadron level where one 
can study the effect of phase transitions on meson and 
diquark properties in mean field approximation, in the 
NJL model applied to quarks, only the order parameters 
of the phase transitions, namely the condensates, can be 
described at mean field level, to study mesons and di¬ 
quarks should go beyond the mean field approximation. 
It is well-known that [27] in the mean field approxima¬ 
tion to quarks together with random phase approxima¬ 
tion (RPA) to mesons in the NJL model, one can obtain 
the hadronic mass spectra and the static properties of 
mesons remarkably well, especially the Goldstone mode 
corresponding to the chiral symmetry spontaneous break¬ 
ing. However, in the pion superfluidity or the color super¬ 
conductivity phase, the structure of the bubble summa¬ 
tion in RPA is much more complicated than that in the 
chiral breaking phase. The correlation among the mesons 
and diquarks becomes very important and leads to the 
Goldstone modes corresponding to the isospin and color 
symmetry spontaneous breaking. In this paper, we study 
analytically and numerically the pion superfluidity, its ef¬ 
fect on the meson properties, and its possible competition 
with the color superconductivity and kaon condensation 
in the frame of standard NJL model at finite temperature 
and baryon, isospin and strangeness chemical potentials. 

The paper is organized as follows. In Section H, we es¬ 
tablish the mean field theory of the NJL model at finite 
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temperature and baryon and isospin densities. In Sec¬ 
tion III, we study the pion superfluidity at zero temper¬ 
ature, zero baryon chemical potential, but finite isospin 
chemical potential, and calculate the chiral and pion con¬ 
densates and thermodynamic functions. In Section IV, 
we study the temperature behavior of the pion super¬ 
fluidity at zero baryon chemical potential and show the 
phase diagrams in the temperature and isospin chemical 
potential plane. In Section V, we investigate the meson 
modes in both normal phase and superfluidity phase in 
RPA. In Section VI, the bosonized version of the SU (2) 
NJL model is discussed at finite isospin chemical poten¬ 
tial with pion condensation. In Section VII, we study 
the baryon chemical potential effect and the competi¬ 
tion between pion and diquark condensations. In Section 
VIII, we briefly discuss the flavor SU{3) NJL model at 
finite isospin and strangeness chemical potentials includ¬ 
ing both pion and kaon condensations. We summarize in 
Section IX. 


II. MEAN FIELD THEORY OF NJL MODEL AT 
FINITE fii 

The standard approach for dealing with the thermody¬ 
namics of a variable particles is via the grand canonical 
ensemble. The key quantity for a system with baryon 
number conservation and isospin number conservation is 
the partition function defined by 

Z(r, HI, HB, V) = . (1) 

Here V is the volume of the system, (3 the inverse tem¬ 
perature /3 = 1/T, hb and hi are the baryon and isospin 
chemical potentials, and B and Is the conserved baryon 
number and isospin number operators. If we take only 
quark held ip as elementary field of the system, the op¬ 
erators can be expressed as 

B ^ ^ J d^xipjoip , 

-fa = ^ / d^:>cipjoT^ip, (2) 

The factors l/W and 1/2 reflect the fact that Nc quarks 
make a baryon and quark’s isospin quantum number is 
1/2. In the imaginary time formulism of finite temper¬ 
ature held theory, the partition function can be repre¬ 
sented as 

ZiT,Hi,HB,V)= J[diP][diP]eIo‘^-I‘‘'<^+^l^^<>^) , (3) 
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The flavor SU{2) NJL Lagrangian density is defined 
as [25-29] 


C = ip 


Too) Ip+ G {ipipy -I- {ipi^zTipy 


( 5 ) 


with scalar and pseudoscalar interactions corresponding 
to a and tt excitations, or equivalently, 


C = ip - Too) Ip 


G 


■ T ,\2 


( 6 ) 


(ipip) + {ipijbTsip) + 2 {lpij5T+1p) {lpij5T-1p) 


with 7r_|_ and 7r_ excitations instead of tti and 7r2, and 

T ± = (n ± iT 2 ) / V ^. 

The above Lagrangian density has the symme¬ 
try [/b(1) 0 S'C/7 (2) 0 S'[/^(2) corresponding to baryon 
number gauge symmetry, isospin symmetry and chiral 
symmetry, respectively. However, in the presence of a 
nonzero isospin chemical potential, the isospin symme¬ 
try SUi{2) breaks down to 17/(1) global symmetry which 
is related to Bose-Einstein condensation of charged pions, 
and the chiral symmetry SUa{2) breaks down to Uia{3) 
global symmetry which is associated with the chiral con¬ 
densation of sigma meson. Therefore, in the case with 
Hi 7 ^ 0, chiral symmetry restoration at finite tempera¬ 
ture and chemical potentials means only degenerate of 
a and ttq meson masses. At zero baryon chemical po¬ 
tential, the “Fermi surfaces” of the u{d) and anti-(i(M) 
quarks coincide. Therefore, at sufficiently high hi > ^ 
the condensate of u and anti-d quarks is favored, and at 
sufficiently high hi < ^ the condensate of d and anti-u 
quarks is favored. Since at low isospin chemical potential 
(and low temperature and low baryon chemical potential) 
there is no deconfinement, the condensates can be iden¬ 
tified as Bose-Einstein condensates of charged pions, 7r+ 
and TT~. Introducing the chiral condensate, 

(iPiP) = a = au + ad , 

Uu = (uu) , 

ad = (dd) , (7) 


and two pion condensates, 


{ipij^T+ip) = V2{uij5d) = 7r+ 
{'ipijiiT-'ip) = V2{dij5u) = Tr~ 




( 8 ) 


where C is the Lagrangian density of the system, and /t 
the quark chemical potential matrix in flavor space jl = 
diag(//„, Hd) with the u and d quark chemical potentials, 


a nonzero condensate a means spontaneous chiral sym¬ 
metry breaking, and a nonzero condensate tt means spon¬ 
taneous isospin symmetry breaking. The phase factor 6 
related to the condensates 7r_|_ and 7r_ indicates the di¬ 
rection of the 17/(1) symmetry breaking. If the system 
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is in global thermal equilibrium, 0 is a constant and it 
does not change any physical result. For convenience we 
choose d = 0 in the following. As is well-known and will 
be seen later, there are corresponding Goldstone modes 
in the region of spontaneous chiral breaking and the re¬ 
gion of spontaneous isospin breaking. 

Defining the meson fluctuations 6 a, Sj^g, 67 ^+, 6 t^_ as 

Ip-ip = a + 6a , 

'iphuTS'ip = S^g , 

-Ppi^zT+ip = 7r+ -I- 6^_^_ , 

'ipi^^T-'ip = TT~ + , (9) 


and keeping only the linear terms in 6 , the partition func¬ 
tion in mean held approximation is simplified as 

= J[diP][di;]eIo I , ( 10 ) 

with the mean held Lagrangian density 

^mf = ^ - (mo - 2Gcr) -I- /xyo -b 2GTrij^Ti] ip 

-G{a^+T:^) , (11) 


from which the inverse quark propagator matrix in the 
flavor space as a function of quark momentum can be 
derived directly. 



+ fj,ulo - Mg 
2fG7r75 


2fG7r75 

7^*^ -b ^d7o - Mg 


( 12 ) 


with the effective quark mass 


Mg = -mo — 2Ga . 


(13) 


Using the method of massive energy projectors [34], we 
obtain explicitly the mean field quark propagator 


o /7 \ _ f ^uu{k} Suc[{k) 


with the four matrix elements 

fco + -S'fc + d-d 


^uu^k) — 


Sddik) = 


Sudik) = 


Sdu{k) = 


(fco 

-EZ){ko + E^) 


ko ~ Ek -b /id 

(fco 

- E+){ko + E+) 


ko - Ek+ fdu 

(fco 

-EZ){ko + E^) 


ko + Ek -b /ill 

(fco 

- U+)(fco + A+) 


2iG'K 

(fco 

-EZ){ko + E-) 


2zG7r 

(fco 

-U+)(fco + A+) 


2zG7r 

(fco 

-EZ){ko + E^) 


2iGn 

(fco 

- E+)iko + E+) 


A-7o 


A+7o 


A-75 


A-H75 


(14) 


(15) 


where are effective quark energies 

_ pit ^ 


Et = ^{Ek ± 1^7/2)" + 4G%2 , 
Ufc = ^|kp+M2 , 
and A± are the energy projectors 

A,(k) = 1 (1 ± 


( 16 ) 


(17) 


A. Gap Equations 


The quark propagator (14) is the background of calcu¬ 
lations for quarks in mean field approximation and also 
for mesons in RPA. From the definitions of the chiral and 
pion condensates (7) and (8), it is easy to express them 
in terms of the matrix elements of the quark propagator, 

o-u = -Ng J [f5„u(A:)] , 

CTd = -Ng J [iSdd{k)] , 

7 r = NgJ [{Sudik) + Sduik)) 75 ] , (18) 

where the trace Tr^j is taken in Dirac space and the 
four momentum integration is defined as f = 

iT f (1^ finite temperature. Obviously, in the 
study without considering color superconductivity, the 
color degrees of freedom in the NJL model is trivial, and 
the trace in color space simply contributes a factor JVg. 
Performing the trace in spin space (the trace including 
the energy projectors are presented in Appendix A) and 
the Matsubara frequency summation we have 


\f{EZ) + f{-E^) - f{E^) 


J (27r)3 Ek 

- f{-Et) + (fiEZ) - /(-A+)) 

Ek 

+ {f{Et) - /(-A+)) 

1+7, I Ek [dj /2 


+ fi-Et) + 


e: 


(fiEZ) - /(-A+)) 


(/(P) _ /(-U+)) 


EP 


= -ANgGir 


1 


(27r)3 


Ek 


(19) 
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with the Fermi-Dirac distribution function 


can be evaluated with the effective quark energies 


fix) 


1 

+ 1 ■ 


( 20 ) 


This group of gap equations has a symmetry that it 
is invariant under the transformations (t„ ^ cr^, (7^ ^ 
CTu,/x/ ^ —Hi- Therefore, we can only concentrate on 
the case fij > 0 and the results for the case fij < 0 can 
be obtained easily. It is also necessary to note that for 
fiB = 0 or fj,i = 0, there is always cr„ = (7^, since in this 
case the chemical potential difference between u and u is 
the same as the difference between d and d. 

In any case we can combine the first two gap equations 
into one determining the chiral condensate cr, 


d^k 2NcMq\Ek- 


+ 


J (2^)3 Ek 

Ek + M//2 


E7 


{f{EZ) - fi-E-)) 


Et 


(fiEf) - fi-E+)) 


( 21 ) 


We now consider the limit of /i/ = 0. In this case the 
gap equations are reduced to 


r 


TT 1 + 
= 0 . 




Only in chiral limit, mo = 0, there is possibly nonzero tt, 
and the gap equation 


1 + 


/ 


= 0 


d^k 8N^G 

(27r)3 yfcJ:^G2(a2+7r2) 


{f{EZ)-fi-E-)) 


(23) 


determines the combination of the two condensates cr^ + 
TT^. In real world with nonzero quark mass, mo ^ 0, the 
pion condensate is forced to vanish, tt = 0 at /i/ = 0, and 
a satisfies the well-known NJL gap equation 


a — 


/ 


= 0 . 


d3k ANcMg 
(2^)3 Ek 


fiEk - ^) - fi-Ek 



B. Thermodynamics 


The thermodynamic potential in mean field approxi¬ 
mation 


n = -L\nZ = -^TT\nS-) + G{a^+7T^) (25) 


7 o-d,7r) = +G(a2 + 7r^) 

-IV. f S^T \ez -E-+Ef-E+ 


(26) 


(277)3 L 


-k2T 


(ln( 


1 -I- e 


-EZIT 


).!„( 


l-ke^+ 


E-/T 


+ ln(l + e-^-/^)+ln (l + ) 


The gap equations (19) for the condensates are equiv¬ 
alent to the extremum condition of the thermodynamic 
potential, 


dfl 

da. 


= 0 , 


dn 

dad 


= 0 , 


dn 

diT 


= 0 . 


(27) 


Note that when there exist multi roots of the gap equa¬ 
tions, only the solution which satisfies the minimum con¬ 
dition 


d^n 

dcrl 


> 0 , 


d^n 


> 0 


d‘^n 

977^ 


> 0 


(28) 


is physical. 

Once n is known, the thermodynamic functions that 
measure the bulk properties of matter can be obtained. 
For an infinite system, these are the pressure p, the en¬ 
tropy density s, the charge number densities ub and nj, 
the flavor number densities and Ud, the energy density 
e, and the specific heat c, that are defined as 


p = -n , 
dn 


® dT ’ 


dn 

dn 

44s = a 

dpB 


dn 

dn 

44“ — "q > 
d pu 

44d = - 

dpd 


e = —p + Ts + piui + Pb^b 1 
de 


(29) 


The flavor number densities = (ujou) and Ud = 
{d'jod) can also be obtained directly from the matrix el¬ 
ements of the quark propagator. 


riu = -Nc 

= Nr, I 


f d^k 

(277)4 

(i3k 


Tru [iSuu{k)jo] 


(27r)^ L' 


f{EZ)+fi-E-) + fiEf) 


+ fi-Et) + {f{EZ) - fi-E^)) 

Ek 

- (fiEf) - fi-E+)) - 2 


Et 


nd = -Nc 


d^k 


'Tr£, [iSddik)jo] 
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= Nr 


d^k 


+ fi-K) - 


■[f{EZ)+f{-E-) + fiE^) 
^^^^{f{EZ)-f{-E-)) 


+ {f(E+) - fi-Et)) - 2 


(30) 


Each pure number density is the difference between 
the corresponding quark number density and antiquark 
number density, 


= n+ - n„ , 

nd = n+ -n~ . (31) 


With the help of the positive and negative energy pro¬ 
jectors A+ and A_, we separate from Uu^d, 


= N 


d^k r 


(27r)3 L 
Ek — iiil"2 


El 


f(EZ)Ef{-E-^) 

{fiEZ)-fi-E-)) 


n„ = —N, 


d^k r 


(27r)3 L 
Ek + M//2 


f{El) + /(-E+) 


K 


nj = Nr 


d^k 


{2'kY L 


{f{Et) - f{-EX)) - 2 


f{Et) + /(-E+) 


+ (/(^;^) - f{-EX)) 


n. = -N, 


(i3k r 


(27r)3 L 
Ek — IJ^il‘2. 


f{EZ)+f{-El) 


El 


{f{EZ)- f{-E-))-2 


(32) 


It is easy to obtain the relations between the flavor num¬ 
ber densities and the charge number densities, 

= ^(n„ - nd) , 

nB = ^{nuEnd). (33) 


III. PURE ISOSPIN EFFECT IN MEAN FIELD 
APPROXIMATION 

To solve the gap equations and calculate the thermody¬ 
namic functions numerically, we should first fix the model 
parameters. Because of the contact interaction between 
quarks that is introduced in the model, there is of course 
no confinement. A further consequence of this feature is 
that the model is non-renormalizable, and it is necessary 
to introduce a regulator A that serves as a length scale in 


the problem, and which can be thought of as indicating 
the onset of asymptotic freedom [41], In the following 
we take a hard three-momentum cutoff A, which and the 
other two model parameters, the effective coupling con¬ 
stant G and the current quark mass toq, can be fixed by 
fitting the pion mass m,r = 0.134 GeV, the pion decay 
constant /,r = 0.093 GeV, and the quark condensate den¬ 
sity (Tu = (Jd = (—0.25 GeV)3 in the vacuum. They are 
[31] mo = 0, A = 0.65 GeV and G = 5.01 GeV“^ in the 
chiral limit, and toq = 0.005 GeV, A = 0.653 GeV, and 
G = 4.93 GeV“^ in the real word. 

In this section we concentrate on the pure isospin effect 
with T = /Zb = 0. In this case the difference between the 
u quark and d quark condensates disappears and we need 
to calculate the total chiral condensate a = 2(j„ = 2ad 
only. 


A. Pion and Chiral Condensates 


At zero temperature and baryon chemical potential the 
gap equations for the chiral and pion condensates are 
reduced to 


fj -j- 2NrAdn 


TT 


1 - ANrG 


d^k I 

r d^k / 

W V 


/ Ek — III 12 

V E^ 

El ^ EX)_ 


Ek + /‘2\ 

Et ) 

= 0 . 


= 0 , 

(34) 


It is clear that the order parameter tt = 0 is always a 
solution of the second gap equation, it corresponds to 
the phase with isospin symmetry at low isospin chemical 
potential. The other order parameter cr describing chiral 
properties in this phase is /i/ independent, 


This means that the chiral condensate keeps its vacuum 
value in the isospin symmetric phase with tt = 0. 

At the critical isospin chemical potential /ij where the 
isospin symmetry starts to break spontaneously and the 
pion condensate appears, the solution tt = 0 should sat¬ 
isfies the equation 


1 - mrG 


d^k 1 El 
(27r)3 Ek El - injf /4 


(36) 


From the comparison with the well-known pole equation 
[26-31] determining the pion mass in the vacuum, 


1 - 8NrG 


d^k 1 El 

^^k El-ml!A 


(37) 


where the quark mass Mq hidden in the quark energy 
Ek is the same as that in (36) and controlled by the gap 
equation (35), we find explicitly that the critical isospin 
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chemical potential f/f for the pion condensation phase 
transition at T = /xb = 0 is exactly equal to the pion 
mass in the vacuum, 


. 


(38) 


For /i/ > /Xj the isospin dependence of the two order 
parameters is governed by the first equation of (34) and 


1 - mM 


(^ 




(39) 


The condensates a and tt scaled by the chiral conden¬ 
sate (To in the vacuum are shown in Fig.(l) as functions of 
/X/. In chiral limit, pions are Goldstone particles of chiral 
symmetry spontaneous breaking, = 0. Therefore, at 
/X/ = 0 the pion condensate jumps up from the vacuum 
value TT = 0 to a finite value, and the chiral condensate 
jumps down from the vacuum value (Jq = 2(—0.25 GeV)^ 
to zero and then keeps zero in the region /x/ > 0 due to 
the constraint of the gap equations (34). These sudden 
changes at the critical point /Xj = 0 means that in chiral 
limit, no matter how small the isospin chemical potential 
is, the chiral symmetry is restored and pion Bose-Einstein 
condensate occurs. In the real world with nonzero cur¬ 
rent quark mass, the critical point is at /if = m,r = 0.134 
GeV and the both chiral and pion superfluidity phase 
transitions are of second order. When the isospin chemi¬ 
cal potential is small, the vacuum state is not disturbed. 
Only when it exceeds the pion mass, the pion condensate 
goes up with /x/ and at the same time the chiral conden¬ 
sate drops down with /x/. The two condensates coexist 
in a wide region. Different from the result obtained with 
chiral perturbation theory [8] where the pion condensate 
at finite /x/ is always less than the chiral condensate in the 
vacuum, tt/cto < 1, tt in our calculation can be larger than 
(To at some intermediate isospin chemical potential. The 
pion condensate tends to zero at about /x/ = 1.75 GeV. 
While the value of the endpoint of pion condensation in 
our calculation is closely related to the finite momentum 
cutoff A, the disappearance of pion condensation at suffi¬ 
ciently high isospin chemical potential can be understood 
in the frame of asymptotic freedom of QCD. 

What is the isospin effect on the chiral condensate if 
we do not consider isospin spontaneous breaking? In this 
case, the only condensate, the chiral condensate cr is con¬ 
trolled by the gap equation 


(T -j- 


d^k 1 
( 2^)3 Ek 


&{Ek-\fii\/2)=0, (40) 


the isospin chemical potential dependence of a is similar 
to its well-known baryon chemical potential dependence. 
In chiral limit, the phase transition point is at 




2WA2- 


2GNc 


(41) 



0 0.5 1 1.5 



FIG. 1. The chiral and pion condensates a and tt, scaled 
by the chiral condensate (Tq, as function of isospin chemical 
potential fii in the chiral limit (upper panel) and real world 
(lower panel) at T = hb ~ 0. 


B. Dispersion Relations 

The dispersion relations of quasi-particles are given by 
the poles of the quark propagator (14), 

u;i(|k|) = EZ(|k|) , u;2m =-E-i\k\) , 
u; 3 (|k|)=^;^(|k|) , u; 4 (|k|) =-if+(|k|) . (42) 

In the case with /xb = 0, it is easy to see the symmetry 
of oxi = — 0 X 2 and 0 x 3 = — 0 x 4 . The dispersion relations 
in this case are shown in Fig.(2) at T = 0. In both 
the isospin symmetric phase with /x/ = 0.1 GeV < 
and the symmetry breaking phase with /x/ = 0.4 GeV 
> itItt, any quasi-particle excitation needs energy. In the 
normal phase with tt = 0 , the chiral condensate cr opens 
the gap for the quasi-particle excitations, and in the pion 
superfluidity phase with /x/ > niTr, the gap is mainly due 
to the pion condensate tt yf 0. In the normal phase, 
the minima of the quasi-particle energies are located at 
|k| = 0 , but in the superfluidity phase, the location of the 
minima of 0 x 1 and CO 2 are shifted to a finite momentum 
when fij is large enough. 
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FIG. 2. The dispersion relations of the quasi-particles 
uii, i = 1, 2, 3,4 in the normal phase with fii = 0.1 GeV (up¬ 
per panel) and pion superfluidity phase with jj,i = 0.4 GeV 
(lower panel) at T = hb = 0. 

The quark occupation numbers n+(|k|), n“(|k|), 
nj'(|k|) and n))'(|k|) in momentum space are just the in¬ 
tegrated functions in (32). At T = = 0 they are 

reduced to 

= njj-dkl) = 1 - , 

«J(|k|) = nCdkl) = 1 - . (43) 

In the normal phase we have n+dk|) = n“(|k|) = 
n)[(|k|) = n^j’dkj) = 0, since the ground state in this 
case is just the vacuum state and there are no quark 
and anti-quark excitations. In the superfluidity phase, 
there are excited quarks and anti-quarks, the occupation 
numbers are shown in Fig.(3) as functions of momentum. 
Note that the maximal occupation number is 2 because 
of the spin degenerate. 


C. Bulk Properties 


We first calculate the isospin density. At T = /is = 0 
it is 


m = Nc 


r 

Ek+fJ^i/‘2 Ek — Hi/2 

J (27r)3 

y Et E^ \ 


(44) 


the numerical result, scaled by the normal nuclear den¬ 
sity no = 0.17/fm^, is shown in Fig.(4). Again in the 


normal phase with /ij < the unchanged ground state 
leads to a zero isospin density. Only when fj,j exceeds 
the vacuum pion mass, nj increases monotonously with 
fij. This nonzero net isospin density is due to the Bose- 
Einstein condensation of charged pions which results in 
different flavor densities. To see this clearly, we plot the 
flavor densities, scaled by the normal nuclear density no, 
as functions of isospin chemical potential in Fig.(5). In 
the superfluidity phase the relation = n'^ > n'^ = n~ 
leads to the net isospin density, and then the number of 
7r_|_ in the system should be larger than the number of 
7r_. 




FIG. 3. The quark occupation numbers n)))jj(|k|) as func¬ 
tions of momentum in pion superfluidity phase with fii = 0.2 
GeV (upper panel) and jii = 0.4 GeV (lower panel) at 
T = hb = 0. 



FIG. 4. The isospin density m, scaled by the normal nu¬ 
clear density no, as a function of /ij at T = /rs = 0. 
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FIG. 5. The flavor densities scaled by the normal nu¬ 
clear density no, as functions of /ij at T = hb = 0. 


Physically, the thermodynamic potential 
ri(r, hb, \o-u(T, fii, cFd{T, fii, ^b), 7r(r, fii, hb)) 
corresponds to the pressure except for a sign, and only 
the pressure relative to the physical vacuum can 

be measured. The physical vacuum is defined to be 

= f}(0,0,0|a„(0,0,0), ad(0,0,0),^(0,0,0)) , (45) 


which for the mean field approximation corresponds to 



(M,(0,0,0)-mo)2 

4G 


-4iVe 


_fk 

(27r)' 


rEk ■ 


(46) 


We introduce the rescaled thermodynamic potential 
n{T,fXl,fJ,B\<Ju,ftd,Tr) = n{T,fXi,fJ,B\<Ju,(Td,Tr) - > 

(47) 

the measurable pressure p and the energy density e are re¬ 
lated to it. Of course, the redefinition of O will not change 
the number densities and the entropy density since they 
are the derivatives of 14 with respect to p, and T. 

The pressure and energy density are shown in Fig.(6) 
at T = pb =0. From the redefinition (47) they are 
zero in the vacuum and keep zero in the normal phase, 
since the isospin effect is not strong enough to disturb 
the vacuum state. In the pion superfluidity phase, the 
ratio of p to e which describes the equation of state of 
the system goes up with increasing pi first and then gets 
saturated, p/e ~ 0.7, at pi ~ 0.5 GeV. This behavior 
indicates that the system is far from an ideal gas which 
satisfies the equation of state p = |e. 

To conclude this subsection, we remark that the isospin 
number density n/ is perhaps a more physical variable 
than the isospin chemical potential pi, since nj can be 
directly measured. The chiral and pion condensates cr 
and TT, scaled by the chiral condensate (Tq in the vacuum, 
the pressure p and energy density e, and the ratio of p to e 
are shown again in Figs.(7) and (8). Different from their 
///-dependence, Figs.(l) and (6), where the isospin sym¬ 
metric and symmetry breaking phases are shown simul¬ 
taneously and are separated from each other at pi = m^, 


the isospin symmetric phase is now contracted at a point 
nj = 0. The isospin density corresponding to the satu¬ 
rated equation of state is about 5 times the normal nu¬ 
clear density. 




and the ratio of p to e (lower panel) as functions of pi at 
T — pB = 0. 



FIG. 7. The chiral and pion condensates a and tt, scaled 
by the chiral condensate (tq in the vacuum, as functions of 
isospin density ni at T = pB =0. 










FIG. 8. The pressure p and energy density e (upper 
panel) and the ratio of p to £ (lower panel) as functions of 
isospin density n/ at T = /is = 0. 

IV. TEMPERATURE EFFECT IN MEAN FIELD 
APPROXIMATION 

We now study the temperature behavior of the chiral 
and pion condensates and the thermodynamic functions, 
and discuss the phase diagram in the T — fxj plane in 
mean field approximation at ps =0. 

A. Chiral and Pion Condensates 

The temperature dependence of the two condensates 
a and tt, again scaled by the chiral condensate ctq in the 
vacuum, at a fixed isospin chemical potential pi = 0.05 
GeV is shown in Fig.(9) in the chiral limit. In the chiral 
limit any small isospin chemical potential will force the 
broken chiral symmetry to be restored and the isospin 
symmetry to be broken at T = 0. Therefore, the chi¬ 
ral condensate keeps zero, and the pion condensate is 
continuously melted in the hot medium. Finally at a 
critical temperature Tc ~ 0.18 GeV the pion condensate 
vanishes and the isospin symmetry is restored. At small 
isospin chemical potential the value of the critical tem¬ 
perature for isospin restoration is almost the same for 
chiral restoration in the study without considering pion 
condensation [27,31]. In Fig. (10) we show the pion and 
chiral condensates as functions of T at pi = 0.15 GeV 


and 0.2 GeV. For /t/ > tUtt, the pion condensate is al¬ 
ready nonzero at the beginning, then drops down due to 
the temperature effect, and finally disappears at a criti¬ 
cal temperature Tc- Like the standard BCS theory, the 
critical temperature and the pion condensate at T = 0 
satisfies the linear relation, 

rc(Ai7)~ 0.577r(r=0,/r/) . (48) 

The relative strength of the two condensates at the be¬ 
ginning depends on the isospin chemical potential, a{T = 
0) > 7r(r = 0) at small and tt{T = 0) > a{T = 0) at 
large /i/. Qualitatively different from the study on chi¬ 
ral symmetry restoration without considering pion con¬ 
densation where the chiral condensate decreases with 
increasing temperature monotonously, here a goes up 
continuously in the coexistence region of the two con¬ 
densates, and drops down only in the region where the 
isospin symmetry is restored. However, we found that 
when Pi is close to the critical value pj = say 
< Pi < 0.2 GeV, the total condensate 
(the dashed line in Fig.(10), behaviors in a similar way 
like the a in the normal study without considering pion 
condensation. 



0 0.05 0.1 0.15 0.2 

T (GeV) 

FIG. 9. The chiral and pion condensates a and tt, scaled 
by the chiral condensate (tq in the vacuum, as functions of 
temperature T at ps = 0 and pi = 0.05 GeV in the chiral 
limit. 

B. Number Densities 

The isospin density ni and flavor densities scaled 
by the normal nuclear density ng, are shown in Fig.(11) 
as functions of isospin chemical potential pi at pB = 0 
but T = 0.1 GeV. At finite temperature, there is thermal 
excitation of particles. Since T = 0.1 GeV is not very 
high, the critical isospin chemical potential /tj for pion 
condensation is still approximately the vacuum pion mass 
niTr- In the region of pi < pf, the small but finite isospin 
density is purely due to the thermal excitation, while in 
the region of pi > pj the rapidly increasing density is 
dominated by the Bose-Einstein condensation of charged 
pions. The thermal excitation depends strongly on the 
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isospin chemical potential carried by the particle. At 
fj-i < the densities of quarks with positive isospin 
chemical potential increase with fxj, while the densities of 
quarks with negative isospin chemical potential decrease 
with /i/. 




by the chiral condensate (tq in the vacuum, as functions of 
temperature T at hb ~ 0 and p/ = 0.15 GeV (upper panel) 
and Hi = 0.2 GeV (lower panel) in the real world. The dashed 
lines are for the total condensate (a^ + /o-q. 

C. Phase Diagram in T — hi Plane 

In the chiral limit when mo = 0, the solution of the last 
gap equation of (19) for tt with cr„ = (Td = 0 separates the 
region of isospin symmetry breaking with tt 7 ^ 0 from the 
region of symmetry restoration with tt = 0. The phase 
transition line demarcating these two regions is given in 
the upper panel of Fig. (12) in the T—hi plane for hb = 0. 
Since an infinite small isospin effect can restore the chiral 
symmetry, there is cr = 0 in the whole plane in the chiral 
limit. In the real world with nonzero current quark mass, 
the system is in isospin symmetric phase not only at high 
temperature and/or high isospin chemical potential but 
also at low isospin chemical potential with hi ^ 

Since the chiral symmetry phase transition is not well 
defined in the case with nonzero current quark mass, we 
did not show the chiral phase transition line. In principle, 
there is cr 7 ^ 0 in the whole T — hi plane. 



Hi(GeV) 



FIG. 11. The isospin density ni (upper panel) and flavor 
densities ^ (lower panel), scaled by the normal nuclear den¬ 
sity no, as functions of isospin chemical potential hi st hb = 0 
and T = 0.1 GeV. 




FIG. 12. The phase diagram in the T — hi plane in the chi¬ 
ral limit (upper panel) and real world (lower panel) at Hb = 0. 
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V. MESONS AT FINITE ISOSPIN DENSITY 

We now investigate meson properties at finite isospin 
chemical potential. In the NJL model, the meson modes 
are regarded as quantum fluctuations above the mean 
field. The meson modes can be calculated in the frame 
of RPA [27]. When the mean held quark propagator is 
diagonal, for instance the case with only chiral condensa¬ 
tion, the summation of bubbles in RPA selects its specihc 
channel by choosing at each stage the same proper polar¬ 
ization function, a meson mode which is determined by 
the pole of the corresponding meson propagator is related 
to its own polarization function IiMM{k) [30,31] only, 

l-2GnMM(fc) =0 . (49) 


then directly use the method developed in [30,31] to cal¬ 
culate the masses of ct, 7r+, 7r_ and ttq. In the superfluidity 
phase with nonzero tt and off diagonal matrix elements 
in the mean held quark propagator, the collective exci¬ 
tations of the system will in principle not be the original 
meson modes in the vacuum but their mixture. Espe¬ 
cially it is expected that there will be a Goldstone mode 
in the superhuidity phase corresponding to the sponta¬ 
neously isospin symmetry breaking. 

In the following we will study the meson properties in 
the normal phase(7r = 0) and superhuidity phase(7r ^ 0) 
separately. 

A. Meson Properties in Normal Phase 


However, for the quark propagator with oh-diagonal el¬ 
ements, like the cases of rj and rj' meson spectrum [27], 
pion superhuidity considered here, and color supercon¬ 
ductivity in Section VII, we must consider all possible 
channels in the bubble summation in RPA. Using matrix 
notation, the ehective quark propagator (14) leads to the 
meson modes determined by 

det(l - 2Gn(A:)) = 0 , (50) 

with the polarization function matrix 


In the vacuum, due to isospin symmetry, one can re¬ 
gard TTi, tt2 or 7r_|_, 7r_ or any other linear combination of 
TTi, tt 2 as eigen collective modes. In the normal phase with 
vanished pion condensate, the polarization function ma¬ 
trix is diagonal with respect to the basis (ct, 7 r_|_, 7 r_, 713 ) 
but not diagonal with respect to the basis (cr, tti , 712 , TTa). 
This means that in the normal phase, the eigen collective 
modes are cr, 7r+, 7r_, ttq. Thus the dispersion relation for 
each mode is determined through the equation 

l-2GnMM(A:) =0 . (54) 


1 -2Gn = 

/1 - 2Gn,, 

-2Gn,,+ 

-2Gn,,^_ 

-2Gn,,„ 


1 - 2 Gn,+,^ 

-2Gn,+,_ 

— 2Gn.n._^,ro 

-2Gn^_,, 

-2Gn,_,+ 

1 - 2 Gn,_,_ 

— 2Gn,r_7ro 

\ ~2Gn,ro(T 

— 2Gn,ro'7r+ 


1 — 2Gn.n.g,ro 


Here the polarization functions defined as 


H 


MM 


(fc) = 


d'^p 


Tr {r^Srufip + k)TM'Smfip)) 


(52) 


with the vertexes 


Tm = 


r 


* _ 

M — 


1 M = a 

iT+j 5 M = 7r+ 

zt_75 M = tt- 

*T375 M = TTq 

1 M = a 
zt_ 75 M = 7r_|_ 

zt+ 75 M = 7r_ 
iT375 M = ttq 


(53) 


are shown in Appendix B. Here the trace Tr is taken in 
color, flavor and Dirac spaces. 

As will be seen in the following, in the normal phase 
with vanished pion condensate, the polarization matrix is 
diagonal at finite isospin chemical potential and one can 


( 5 l)rhe meson mass is defined as the root fcg at k = 0. 

\ Doing the trace in flavor and color spaces, we obtain 
the polarization functions in terms of the matrix elements 
of the quark propagator, 

^cr<y{k) = iNc J -^^Ttd(Suu{p + k)Suu{p) (55) 
+ Sddijp + k)Sddip)^ ) 

nTToTTo)^) ~ fA^C j' ^2^)4 {p T k)^^Suu{p) 

+ JbSddip + k)"f5Sdd{p)^ , 

n^+7r+(fc) = -SzA^c J ■^^T'T^D(j5Suuip+k)-f5Sddip)j , 

( 27 r ^4 Trj (^^^Sddip + k)^zSuu{p)^ , 

now the trace Tr^i is taken only in Dirac space. At zero 
temperature and zero baryon density, the polarization 
functions at k = 0 are simplihed as 


n,Ako) = m,J 




H 


d^p 1 

(27r3) Ej, El - PJA ’ 
d^p I El 

772 h , 
0 / 


(fco) > 

.(fco) =4Afc J 


d^p 1 


El 


(27r3) Ep El - {ko + Az/) V4 ’ 
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n^_7r_(^o) = 4A^c J 


d^p 1 




(27r3) Ep El - (fco - M/)74 


• (56) 


From the comparison with their expressions in the vac¬ 
uum with /X/ = 0, we obtain the fij dependence of the 
meson masses, 




Maiiii) = ruo 
= 'm-n 
=m-„ - , 

=771^,+ . 


(57) 


The mesons which carry zero isospin charge keep their 
vacuum masses, while the mass of the meson which 
carries positive (negative) isospin charge decreases (in¬ 
creases) with Hi linearly. These relations hold before 
Hi = TO,r where the pion condensation starts. 

After performing Matsubara frequency summation the 
meson polarization functions at k = 0 are explicitly ex¬ 
pressed as 


d^p 1 El — 


n.4A:o)- 2iV=/(2^3)£:^£:2_fc2/4 

X (/(AI) + f{Et) - fi-E^) - f{-E+)) , 


1 




(2^3) Ep El - kl/A 


n 


ro7ro(7) = -‘2Nc J 

X (/(AI) + f{E+) - f{-E-) - f{-E+)) , 

iko) = m.J 


Tr+TT+v-uy - --c / 

2Ep — Hi E ko 


X 


kl - 4{Ep - hi/2Y 

2Ep + Hi ~ ko 


{f{EZ) - /(-A+)) 


kl - A{Ep + Hi/2y 


{f{E+) - fi-E+)) 




,_{ko)=^N,j 


d^p 

(^ 

2Ep Hi ko 


kl - i^{Ep - hi/2Y 
2Ert E IIT E k\ 


ifiEZ) - fi-EZ)) 


kl - d{Ep E hi/2Y 


if (El) - fi-Et)) 


(58) 


From the comparison of n 7 r+ 7 r+ ai hb = 0 and fco = 0 
with the last gap equation of (19) for the pion condensate 
TT , the phase transition line from tt = 0 to tt 0 in the 
T — Hi plane at hb =0 can be calculated by the equation 


1 - 2Gn 


7r_|_7r_|_ 


(7 = 0 ) = 0 . 


(59) 


Thus the mass is always zero on the phase boundary. 


B. Meson Properties in Superfluidity Phase 

The polarization functions IImmK^) at k = 0 are 
calculated in Appendix B as explicit functions of tem¬ 
perature T, baryon chemical potential hb, and isospin 
chemical potential hi ia general case with nonzero chiral 
and pion condensates. While ct, 7r_|_ and 7r_ are no longer 
collective excitation modes due to the pion condensation, 
there is still no mixing between ttq and other mesons, 

n,roo-(fc) = n^oy+(fc) = nyroyi--(fc) = 0 , (60) 

its mass in the superfluidity phase is determined by 

1 - 2Gn,„,„ (fco = , k = 0) = 0 . (61) 

Taking the polarization function IlTroTro evaluated in Ap¬ 
pendix B and making comparison with the gap equation 
for TT at T = /xb = 0, we get analytically 

M^voidi) = di , Hi > bItt ■ (62) 

The other meson masses are now determined by the 


equation 




( 1 - 2Gn,, 

-2Gn,,^ 

-2Gn,^_ 

det 1 

-2Gn,+, 

I - 2Gn,+,^ 

-2Gn,+,_ 

= 0 

-2Gn,_, 

-2Gn,_,+ 

I - 2Gn,_, 

at k 

= 0, namely. 




ko=M 

(63) 


1 - 2Gn,+,^(fco))(i - 2Gn,_,_ (fco)) - dG^n^ 


7r_|_ TT — 


x(l - 2Gn,,(fco)) - 16G3n,,+ (A:o)n,._ (A:o)n,+._(7) 


-4G"n7_(A:o)(l-2Gn 
x(l-2Gn,_,_(fco)) 


.+.+ (fco))-4G"n 
= 0 . 


2tt2 

(T7r+ 


(7) 


J ko^M 


(64) 


We first discuss the masses in the chiral limit. In the 
chiral limit, the chiral condensate a as well as the effec¬ 
tive quark mass Mq are forced to be zero by any small 
isospin chemical potential. From the explicit form of the 
polarization functions evaluated in Appendix B, the mix¬ 
ing between a and 7r± disappears, Tlcnv± = 0, then a 
is still an eigen collective mode. Since the polarization 
functions for a and ttq are exactly the same in this case, 
Hfja = the mesons tt® and a have the same mass 

in the superfluidity phase, 

McriT,Hi,HB) = M^q{T,Hi,HB) , (65) 

for any T,hi and hb, which reflects the restoration of 
chiral symmetry. In the case of T = /xb = 0 we have 

M„{hi) = M^q{Hi) = Hi ■ (66) 

The other two eigen modes are determined by 
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(1 - 2Gn,^,^(fco))(l - 2Gn,_,_ (fco)) 


-4G^n^ , 

7r_|_7i 


.(ko) 


ko=M 


= 0 


(67) 


From the polarization functions shown in Appendix B, it 
is easy to prove that M = 0 is a solution of (67) corre¬ 
sponding to the spontaneously isospin symmetry break¬ 
ing. We call the massless meson mode as tt^ and the 
other heavy mode as tth- The meson mass spectra in the 
chiral limit are shown in Fig. (13) as functions of p./ at 

T = fj,B = 0- 



0 0.1 0.2 0.3 0.4 

(GeV) 

FIG. 13. The meson mass speetra in the ehiral limit as 
functions of pi at T = pB ~ 0. 


In the real world with nonzero current quark mass, 
there is not only the mixing between the charged pi- 
ons, but also a mixing between cr and the charged pions. 
From the explicit expression of the polarization functions 
shown in Appendix B, we have 

HcTTri = n,r±<T OC V2Mg\2GTT\ , 

n,r+ 7 r_ = OC 4G^7r^. (68) 


To quantitatively see the mixing degree among cr, 7r_|_ and 
7r_, we introduce two dimensionless parameters. 


V2Mq\2GTT\ 

4GV2 

_ dG^TT^ 

“ 4G^ ’ 


(69) 


describing the strengthes of cr — tt mixing and tt — tt mix¬ 
ing. Their isospin dependence is shown in Fig. (14) at 
T = pB = 0. The a — TT mixing is much weaker than the 
TT—TT mixing, since the former is proportional to both chi¬ 
ral condensate and pion condensate, but the latter to the 
square of pion condensate, and in the superfluidity phase 
the chiral symmetry is almost restored. The a — tt mix¬ 
ing is important only around the phase transition point 
pf = m,r) but the tt — tt mixing is strong in a wide region 
after the phase transition. 

The mass spectra calculated through solving the 
Eq.(64) numerically are shown in Fig.(15), together with 
M,ro) Eq.(62) in the superfluidity phase and the masses 


(57) in the normal phase, as functions of pi at T = pB = 
0. Again we obtain the Goldstone mode reflecting cor¬ 
rectly the spontaneous isospin symmetry breaking in the 
superfluidity phase. In the normal phase cr, 7r_|_, 7r_ and 
ttq themselves are eigen modes of the collective excita¬ 
tion of the system, but in the superfluidity phase, except 
for ttq which is still an eigen mode, cr, 7r_|_ and 7r_ are 
no longer eigen modes. Since the masses determined by 
(64) are continuously connected with the masses of cr, 7r_|_ 
and 7r_ at the phase transition point pi = m,r, we call 
them S,n+ and n_, shown in Fig.(15). When pi is high 
enough, the pion condensation will disappear due to the 
asymptotic freedom of QCD, this turning point in the 
NJL model is about p'j = 1.75 GeV, see Fig.(l). Above 
this turning point the mass spectra are again controlled 
by the polarization functions (56) and the pi dependence 
is similar to (57). 



11, (GeV) 

FIG. 14. The mixing strengthes riair and as functions 
of isospin ehemieal potential pi at T — pB = 0- 

The mass spectra in the superfluidity phase calculated 
in the NJL model here are very similar to the result pre¬ 
dicted in chiral perturbation theory [8,12], 

A7n+ = 0 , 

A7n_ = Pi\Jl + 3rn^/pj , 

= PI . (70) 

However, if we neglect the a — tt mixing in the super¬ 
fluidity phase, cr is still the eigen mode and the mass is 
calculated from the pole equation 

l-2Gn,,(fco=M„k = 0) = 0, (71) 

and the other two eigen modes which are linear combina¬ 
tions of 7r_|_ and 7r_ are determined by the equation (67). 
One of the modes is still the Goldstone mode and the 
other is a heavy one. The numerical results for a and the 
heavy one are also shown in Fig.(15) as dashed lines. In 
this case, the broken CAi(l) symmetry will be restored at 
sufficiently high isospin density, indicating by the same a 
and TTo mass. It is clear that the difference between the 
full calculation and the approximation neglecting a — tt 
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mixing is mainly in the narrow region above the phase 
transition point fj-j = rriT^. 



0 0.1 0.2 0.3 0.4 

11, (GeV) 


FIG. 15. The meson masses as functions of isospin chemical 
potential fij at T — fiB = 0 in the real world. The solid 
lines are the full calculation and the dashed lines from the 
approximation neglecting a — n mixing. 


Using the analytical form of 1122 evaluated in Appendix 
B, we find that the above equation is exactly the same 
as the gap equation for the pion condensate. There¬ 
fore, in superfluidity phase with tt 0, there is always a 
Goldstone mode. Note that this massless mode holds in 
the whole superfluidity phase at finite temperature and 
isospin and baryon chemical potentials. 

The Goldstone mode will extremely change the ther¬ 
modynamics of the system. In mean field approximation, 
only quarks contribute to the thermodynamic potential, 
see (26). Since quarks are heavy at low temperature, the 
thermodynamic functions scaled by the temperature, for 
example, p/T'^, e/T"*, s/T^ and c/T^ approach to zero in 
the limit of T ^ 0 [31]. However, when the meson fluc¬ 
tuations are included, the massless meson mode will lead 
to power laws of thermodynamic functions at low tem¬ 
perature, these scaled functions will approach to nonzero 
constants. 


We now analytically prove the Goldstone mode in gen¬ 
eral case in the NJL model. To this end we consider the 
polarization function matrix 

1 - 2Gn = (72) 

l-2Gnoo -2Gnoi -2Gno2 -2Gno3 \ 

-2Gnio 1 - 2 Gnii -2Gni2 -2Gni3 

-2Gn20 -2Gn2i 1 - 2Gn22 -2Gn23 

-2Gn30 -2Gn3i -2Gn32 1 - 2Gn33 

with respect to the basis (cr, tti, 7r2 ,tts). 

From Appendix B, 7r3 does not mix with other modes, 

n 3 n(fc) = n„ 3 (fc) = 0 , n = 0,1,2, (73) 

its mass is determined by its own polarization function 
at k = 0, 


1 - 2Gn33(Ao = = 0 . (74) 

The masses of the other modes are solved from the 
equation 


/ 1 - 2Gnoo 

-2Gnoi 

-2Gno2 \ 


det -2Gnio 

1 - 2Gnii 

-2Gni2 


V -2Gn2o 

-2Gn2i 

1 - 2Gn22 ) 

ko^M 

= 0 



(75) 

at k = 0. The existence of the Goldstone mode means 


that M = 0 is a solution of the above equation. This is 
really true. In fact, setting fco = k = 0 we have from the 
Appendix B 

1102 ( 0 ,0) = n 2 o(o, 0) = ni 2 (o, o) = n 2 i(o, o) = o , (76) 
the Goldstone mode can be checked by the equation 

1 - 2Gn22(fco = 0, k = 0) = 0 . (77) 


VI. PION SUPERFLUIDITY IN LINEAR SIGMA 
MODEL 

In this section we study the linear sigma model at finite 
isospin chemical potential and compare it with the NJL 
model. Since the linear sigma model can be regarded as 
a bosonized version of the NJL model, we want to know 
whether the former can reproduce the results we have 
worked out with the latter with quarks as elementary 
degrees of freedom. It is well-known that one can use 
the bosonization technique to integrate out the quark 
degrees of freedom and obtain the effective Lagrangian 
for the meson bound states in the NJL model [27,43]. In 
the real world the Lagrangian is just the one of the linear 
sigma model with an explicit chiral breaking term, 

^ -f (9^7f)2] -p 1 {2glfl - ml) + tt^) 

- , (78) 

where = g^^qq is the pion-quark-antiquark coupling 
constant, the pion decay constant, and the pion 
mass in the vacuum, g^^ and can be calculated in the 
NJL model in RPA. 

A. Zero Temperature 

Introducing the isospin chemical potential pi corre¬ 
sponding to the third component of isospin charge 

h = J d^yi{'7ridtn2 - n2dt7Ti) , (79) 

and the chiral and pion condensates 
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(o') = C . 

(7r+) = (tt") = , 


V2 


(80) 


we obtain the effective Lagrangian at tree level, 

^eff = -U{^,p) + i (5^cr)^ + {df,TT3f + {dtTTl - ^1^2^ 
+ (9t7r2 + /X/TTl)^ - (Vtti)^ - (V7r2)^ 


Qf2 I ^2 f2 . "‘"IT 2 

•5? + P - + TTT 


25^ 


+ {e+5p^-f^+ 


25? 


+ U" + - /? + 


25^ 


(tt? +7r?) +4^pcr7ri 


- 25 ? (Co+ p7ri) (cr^+7r2) - ^ (ct^+ 71^) , (81) 

with the classical potential 


U((.p) = -i ( 23 lfi - mi) + I (1^ + 


c 2 t- 1 2 2 


(82) 


where the two condensates are determined by minimizing 
the potential, 




(83) 


We have chosen pion condensate to be real as in the NJL 
model. In this section we use C and p standing for the 
chiral and pion condensates in order to avoid confusion 
with the definitions in the NJL model. 

In the vacuum the constraints (83) on the potential 
give ^ = fir and p = 0, and the sigma mass can be read 
out from the quadratic term in the sigma held in the 
effective Lagrangian (81), to? = dp?/? + to?. 

At hnite isospin density the two condensates in the 
ground state satisfy / = 0,p = /^a/1 + 2p?/to? in the 
chiral limit. This is consistent with the NJL model: Any 
small isospin density can force the chiral symmetry to be 
restored in the chiral limit. In the real world, we obtain 
C = /ti-) P = 0 for pi < m.n which is the same as in the 
vacuum, and [5] 


c 


/tt 2 > 

P/ 


P = Itt 



P? / TO? - TO? ’ 


(84) 


for Pi > TO,r- This result is qualitatively consistent 
with the NJL model and the chiral perturbation theory. 
Especially the critical behaviors of the two condensates 


around the phase transition point p‘f = to^^. in the three 
effective models are almost the same. 

The thermodynamic functions can be evaluated from 
the potential U. The isospin density n/, pressure p and 
energy density e are analytically written as 


ni = 


dU 

dpi 


= PIP 


= pjflU-^.+2 ^\ 

V P/ "T-i- 

P = -ui^,p) + uiu,o) 

= (/ + 


2 

€ = -p 


Pi 


■ pini 


1 f2(.2 2. ( P? + 3to? , 3/x2+to? 


• (85) 



^ H|(GeV) 


FIG. 16. The ratio p/e (upper panel) and the chiral and 
pion condensates (lower panel) as functions of isospin chem¬ 
ical potential p/ at T = ps = 0 in the NJL model (NJL), 
linear sigma model (LSM), and chiral perturbation theory 
(CPT). 

If we take the limit rria ^ oo, we reproduce the result 
in the chiral perturbation theory [8,12]. However, the 
finite sigma mass leads to a big difference at large pj. 
At large pi, the ratio p/t approaches I in the chiral per¬ 
turbation theory, but only about 1/3 in the linear sigma 
model and 0.7 in the NJL model, as shown in Fig.(16). 
While the three effective models are very different in high 


15 








density region, they behavior almost the same in the re¬ 
gion close to the critical point, see also the chiral and 
pion condensates shown in the lower panel of (16). 

The meson mass spectrum is easily obtained by diag¬ 
onalizing the quadratic part of the effective Lagrangian 
(81). It is clear that ttq itself is an eigen mode of the sys¬ 
tem, while 7r_|_ and 7r_ are not, like the case in the NJL 
model. A simple algebra calculation gives the result 

M„{ni) = ma- , 

=m^ + , 

=m^- fii , 

( 86 ) 

for fii < and 






2 

7r«^ TT 


1 

(m/) = 


M 

2 






\ 




Mu+ipi) = 0 , 
M^oiPi) = Pi ) 




4ep2 , 


(87) 


for /i/ > ttItt. The numerical result is shown in Fig(17). 
The a — TT mixing here plays also an important role for 
the meson mass spectra. The calculation without consid¬ 
ering a — TT mixing is indicated by dashed lines. We find 
that the mass spectra are qualitatively consistent with 
the NJL model and the chiral perturbation theory. 



FIG. 17. The meson masses as functions of isospin chemical 
potential /i/ at T = 0 in the linear sigma model. The solid 
lines are the full calculation and the dashed lines from the 
approximation neglecting cr — tt mixing. 


B. Finite Temperature 

There are two approaches to treat thermal excitation 
in the linear sigma model, the Hartree-Fock approxima¬ 


tion and the large N expansion. It is well-known that 
a disadvantage of the Hartree-Fock approach in the dis¬ 
cussion of chiral symmetry restoration is the lack of the 
Goldstone mode in the symmetry breaking phase [44]. 
Here we calculate the mass spectra and phase diagram 
in the two approaches at finite temperature and isospin 
density and compare the results with the NJL calcula¬ 
tion. 

Since the sigma model (78) can be regarded as a 0(4) 
model, we adopt the large N expansion method in the 
0{N) sigma model with boson chemical potential. With 
the same method and technics given by Harber and Wel¬ 
don [45], the mass parameter M and condensates ^ and 
p are determined by the following gap equations 

= 2gl [^2 +p^-fl + 2 J'(r, pj, M)] + ml , 

p{M^ -pj)=0, 

- Uml = 0 , (88) 


the thermal excitation is introduced by the function 
d^k 1 


J'(T,/x/,M) = i J 


(27r)3 ^/WTW 


X [2/fc(T, 0, M) + /fc(T, PI, M) + fbiT, -pi, M)] , (89) 


with the Bose-Einstein distribution function 
fb(T,pi,M) = 

It is easy to derive the solution 


(90) 


M{T,pi) = Pi , 

aT,Pi) = U^ , 

Pi 

p{T, pi) = ^p^{Q,pi)-2J'{T,pi,pi) , (91) 


in the pion superfluidity phase, and the coupled equations 




= 2gl[e-fl + 2J'{T,pi,M)] 




(92) 


for M and ^ in the normal phase with p = 0. From the 
above two groups of solutions we can obtain the phase 
diagram in the T — pi plane, shown in Fig.(18). 

The Hartree-Fock approach in the superfluidity phase 
is rather complicated, and there is no Goldstone mode 
at finite temperature in this approach. Here we will not 
present the treatment in the superfluidity phase with p ^ 
0. However, to calculate the phase transition line of pion 
superfluidity in the T — pi plane is much easier. We need 
only the formulas in the normal phase with p = 0. From 
the detailed derivation given in Appendix G, the effective 
meson masses and the chiral condensate in this phase are 
governed by the coupled equations. 
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Ml = 2gl - fl + 3{aa) + 2(Tnr) + (ttotto)) + ml , 
Ml^ = 2gl (^2 - fl + (ctct) + 2(7r7r) + 3(7ro7ro)) + ml , 
Ml^ = 2gl - fl + (aa) + 4(7r7r) + (ttotto)) + ml , 
‘^gli - fl + 3(o-cr) + 2(7r7r) + (ttotto)) + ^ml = f^ml . 

(93) 

To illustrate how to determine the critical temperature 
from the quantities obtained in the normal phase, we con¬ 
sider the isospin density derived from the thermodynamic 
potential shown in Appendix C, 

do. f 

(94) 

Note that this is only the thermal contribution to the 
isospin density and true only in the normal phase. In the 
superfluidity phase we should include the contribution 
from the condensation of charged pions. From the defini¬ 
tion of Bose-Einstein condensation, the phase transition 
point where the occupation number density is divergent 
at zero momentum is determined by 

M^{Tc,gi) = gi ■ (95) 

We have checked that this equation can also be derived 
from the the superfluidity phase in the limit of p = 0. 




T — Hi plane in the chiral limit (upper panel) and the real 
world (lower panel) in the linear sigma model. The solid and 
dashed lines are, respectively, the results of the large N ex¬ 
pansion and the Hartree-Fock approximation. 


The phase diagram in the T — fij plane calculated with 
the large N expansion method and Hartree-Fock approx¬ 
imation is shown in Fig.(18) in the chiral limit (upper 
panel) and the real world (lower panel). It is clear that 
the result of the large N expansion is qualitatively consis¬ 
tent with the NJL model calculation. Especially, in the 
chiral limit, when gj 0 the critical temperature in the 
large N approach is very close to the one obtained in the 
NJL model, while in the Hartree-Fock approximation the 
critical temperature approaches zero, which is certainly 
wrong. 

VII. FINITE BARYON DENSITY AND 
DIQUARK CONDENSATION 

We have investigated the NJL model at finite isospin 
chemical potential gj and temperature T and made com¬ 
parison with the other effective models. We now turn 
to the discussion at finite baryon chemical potential gs- 
It is well-known that [7,27,47] the baryon density effect 
on deconfinement and chiral restoration is qualitatively 
different from the temperature effect: the phase transi¬ 
tion is of first order at high baryon density but of second 
order or even a smooth crossover at high temperature. 
The physics in high baryon density region is also very 
different from that in high isospin density region: the 
spontaneously broken symmetry is of color symmetry at 
sufficiently high baryon chemical potential but isospin 
symmetry at sufficiently high isospin chemical potential, 
the former is described by the diquark condensate (tpip) 
and the latter by the pion condensate tt. We will first 
consider the baryon density effect on the pion conden¬ 
sation without considering the color superconductivity, 
and then discuss the competition between the pion su¬ 
perfluidity and color superconductivity in the gs — gi 
plane. 

A. Finite Baryon Density Effect without Diquark 
Condensate 

As we indicated in the Section B, when the two chem¬ 
ical potentials gs and gj are both nonzero, the u and 
d quark condensates are different from each other, see 
(19). It is clear that the isospin chemical potential gi 
should be large enough to guarantee the phase transition 
of pion superfluidity at finite baryon chemical potential 
gs and temperature T, at least it should be larger than 
the critical value m^^ at T = gs = 0. In our numerical 
calculation we took it as 0.2 GeV. The gs dependence 
of the chiral and pion condensates, again scaled by the 
chiral condensate (Tq in the vacuum, at fixed T and gi 
is shown in Fig.(19). At zero temperature, the u and d 
quark condensates cr„ and Ud are almost the same, and 
all the three condensates, (j„, Od and tt keep their vacuum 
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values in the pion superfluidity phase. At a critical value 
fiB ^ 0.85 GeV, a first order phase transition happens, 
the pion condensate jumps down from its vacuum value 
to zero, and the u and d quark condensates jump up from 
their vacuum value to some larger values. In the normal 
phase with tt = 0, the two chiral condensates decrease 
with increasing baryon chemical potential, and the dif¬ 
ference between them becomes remarkable. The temper¬ 
ature effect changes the phase transition from first order 
to second order, see the lower panel for T = 0.1 GeV, 
and reduces the critical baryon chemical potential from 
about 0.85 GeV to about 0.65 GeV. Also the tempera¬ 
ture effect leads to remarkable difference between the two 
chiral condensates not only in the normal phase but also 
in the pion superfluidity phase. 




FIG. 19. The u and d quark condensates and pion conden¬ 
sate, scaled by the chiral condensate uo in the vacuum, as 
functions of /is at /i/ = 0.2 GeV and T = 0 (upper panel) 
and T = 0.1 GeV (lower panel). 

The temperature dependence of the three condensates 
is indicated in Fig. (20) at fixed chemical potentials. Sim¬ 
ilar to the case [30,31] of investigating chiral symmetry 
restoration without pion condensation, the temperature 
effect here results in a second order phase transition of 
pion superfluidity. The phase diagram in the T — hb 
plane at fixed fxj is shown in Fig. (21), which is very sim¬ 
ilar to the chiral phase transition line at /i/ = 0 [31]. 
The phase transition is of second order in high temper¬ 
ature region and of first order in high baryon chemical 


potential region, and the tricritical point which connects 
the first and second order phase transitions is located at 
T = 0.045 GeV and ^lB = 0.78 GeV. 



sate, scaled by the chiral condensate (Tq in the vacuum, as 
functions of T at /ij = 0.2 GeV and /is = 0.6 GeV. 
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FIG. 21. The phase diagram of pion superfluidity in the 
T — hb plane at jj,i — 0.2 GeV. The solid part of the line 
indicates a first order phase transition while the dashed part 
indicates one of second order. The symbol x is the position 
of the tricritical point. 


B. Competition between Pion and Diquark 
Condensates at Ac = 2 

To study both pion superfluidity and color supercon¬ 
ductivity at finite baryon and isospin chemical poten¬ 
tials, we first consider the two color NJL model with 
both quark-antiquark and diquark channels, 

£ = ■4’ - mo) i^ + Gs + (^* 75 ^'*/')^) 

+ Gd (jp''ij5T2Ti'4;'^ , (96) 

where r| and T 2 are the second Pauli matrix in color and 
flavor spaces, Gs and Gb are the coupling constants in 
the quark-antiquark channel and diquark channel, and 
from the Fierz transformation we have Go = Gs for 

Nc = 2. 
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At sufficiently small m and large /x_b, there should be 
only chiral condensate a and diquark condensate A de¬ 
fined as 


A = (V’°*75'rfr/ \p) . (97) 


Taking the similar way used in Section II, it is straightfor¬ 
ward to derive the quark propagator matrix in Nambu- 
Gorkov space and then obtain the gap equations for the 
two condensates [34,35]. At zero temperature they are 


(7 -t- 


tfk I fEk-fis/Nc 


(2^)3 Ek L E-, 


e(AX-|M/|/2) 


Ek + Mb/A^c 


A 


1 - 
1 

Ea 


©(£;+-Im/I/2) 

d3k 

3 I 17 

A 


= 0, 


&{E^ + \^iI\/2)) 


(98) 


with the energy functions E^ and effective quark mass 
Mq defined as 


Ea = \/(Ek ± + 4G|,A2 , 

Mq = mo- 2Gs<J . (99) 

These gap equations for a and A are exactly the same 
as the gap equations for cr and tt shown in Section II at 
T = 0, if we do the replacement A ^ tt, /is ^ /x/. Mak¬ 
ing comparison with the pion mass equation (37) in the 
vacuum and considering the symmetry between A and 
TT, we obtain the critical baryon chemical potential /x^ of 
color superconductivity and the critical isospin chemical 
potential /xj of pion superfluidity. 


where a, b and c are constants. We can easily obtain a 
first order phase transition line, shown by the dashed line 
in Fig.(22), by minimizing the effective potential. The 
phase diagram calculated here in the frame of NJL model 
is exactly the same as the one obtained with effective 
chiral Lagrangian for Nc = 2 [13]. 




FIG. 22. The phase diagram of pion superfluidity and color 
superconductivity in the — l^i plane at T = 0 in the frame 
of Nc = 2 (upper panel) and A7 = 3 (lower panel) NJL mod¬ 
els. 


Mb = < m^) , 

m5 = , (mb < m ^) ■ 


( 100 ) 


C. Competition between Pion and Diquark 
Condensations at Ac = 3 


Therefore, the normal phase without pion and diquark 
condensations is in the square box /x/ < and mb < 
niTr, shown in Fig.(22). The system is in the phase of 
color superconductivity with tt = 0 and A 0 for /x/ < 
and mb > iti-k and in the phase of pion superfluidity 
with TT ^ 0 and A = 0 for /xb < and m/ > itt-tt- 
In the region of high /x/ and /xb, due to the symmetry 
between the pion and diquark condensates proposed by 
their gap equations, we can conclude that the pion and 
diquark condensates are equivalent in the phase diagram 
in the M/ ~ Mb plane. This can be confirmed by taking 
the bosonized version [42] of the two color NJL model 
with chiral, pion and diquark condensates. The effective 
potential in the bosonized version can be derived as 

V (cr, TT, A) = —a(CT^ -I- TT^ -I- A^) -I- b{a‘^ -I- -I- A^)^ 

- ca - (101) 


In the real world with color degrees of freedom Nc = 3, 
the NJL model including the diquark sector is defined by 
the Lagrangian density 

C = 'tp - mo) -tp + Gs 

+ Gb (viL*756“^'^eyV’i/s) (V’ia*75ec^^e)>|/3) > (102) 

where and totally antisymmetric tensors in 

flavor and color spaces. 

At sufficiently small m/ and large mb, there should be 
only chiral and diquark condensates. Taking the stan¬ 
dard way used in Section II we can derive the quark 
propagator matrix in the Nambu-Gorkov space and then 
obtain the gap equations determining the chiral conden¬ 
sate cr and diquark condensate A defined as 

^ e“^^*75V'i/3) . (103) 
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where it has been regarded that only the first two colors 
participate in the condensate, while the third one does 
not. At zero temperature, the and /i/ dependence of 
the two condensates are governed by the equations 


(7 


1 rEk-HB/i 


(27r)3 Ek L K 


Q{E^ - Im/I/ 2 ) 


A 


+ ^'" p”/^ 9(Et-|wl/2) 

+ i 0 (Ao- - Im/I/ 2 ) + ^e{E+ - |/i,|/ 2 ) 

/ d^k / 1 


= 0 , 


E+ 


eiE+ - \^,j\/2 


= 0 . 


(104) 


Similar to the meson mass equations derived in RPA 
[26,27,31], we can obtain the diquark mass in the vac¬ 
uum by the pole equation 


1 — 8Gd 


d^k / 1 1 \ 

{2Trf \Ek + ttid/^ Ek - mD/2) 

(105) 


At sufficiently small fj-B and large /r/, the system 
should be in the pion superfluidity phase. The chiral 
and pion condensates a and tt are controlled by the gap 
equations shown in Section II. From the comparison of 
them at zero temperature with the pion mass equation 
in the vacuum, and the comparison of the gap equations 
(104) with the diquark mass equation (105), we get the 
conclusion that in the /is — /x/ plane the normal phase 
without pion and diquark condensates is in the rectangle 
defined by 


m KuIt, , 

3 

Ms < -jjriD , (106) 

shown in Fig.(22). Outside this rectangle, we know that 
the system is in color superconductivity phase with A ^ 
0 and TT = 0 for /i/ < m.^. and hb > |tod, and in pion 
superfluidity phase with tt 7^ 0 and A = 0 for ^b < fuzD 
and /i/ > tUtt. As for the phase structure in the region 
with high baryon and isospin chemical potentials, it needs 
to consider the gap equations for the three condensates 
CT, TT and A simultaneously. The dashed line in the lower 
panel of Fig. (22) which separates the phase with tt = 0 
and A 7^ 0 from the phase tt 7^ 0 and A = 0 is just an 
estimation of us. 


VIII. EXTENSION TO SU{i) NJL MODEL 


in the frame of flavor SU{3) NJL model. The Lagrangian 
density is defined as [26-29] 

8 

£ = - TOo)V’ + + ('0*75AaV')^] 

a—0 

- Ar[det 2 /j(l -I- 75 )'!/' -I- det 2 p{l - 75 )V'] , (107) 

where mo = diag{mou,iTiod,'rnos) is the mass matrix of 
current quarks, G and K are coupling constants, and 
the t’Hooft’s determinant includes six-fermion interac¬ 
tion. The three flavor NJL Lagrangian can be brought 
into an effective form similar to (5) of the two flavor NJL 
model by writing the six-fermion interaction in an effec¬ 
tive four-body form [26-29], 

- mo)^ (108) 

8 

+ ^ [G-(V;AV)" + G+(05^75AV)"] 

-b [Go 3 ( 2 /;A° 0 )( 2 /;A^ 2 />) -f Gj3('0i75A°2/>)(0i75A^0)] 

-b -b G;[o(2/;i75A^2/>)(0i75A°0)] 

+ [Go8 (*/’A°0)(2/;A®2/>) -b G[}8(2/;i75A°'0)(0i75A®0)] 

+ [G8 o(*/’A®0 )(0A°2/>) -b Gjo('0i75A®2/>)(V;i75A°0)] 

-b [G;j8(*/’A^0)(*/’A®2/>) -b G;[8('0*75A^0)(0 i75A®0)] 

+ [G83(*/’AV)(*/’A^*/’) + Gj3(2/;i75A®2/>)(0i75A^0)] , 

with the effective couplings 

GJ = G =F 2 A*cAr(cr„ + UdE ) 

Gf =G± =Gf =G±iA,iLa, , 

G^ = Gi = G ± ^N,Kad , 

G± =G± =G±iA,iLa„ , 

G^ = G ± —NcK(2au -b 2ad — (Ts ) , 

n 

Gqs = = ±^A^cA:(cr„ + ad- 2as) , 

^38 “ ^83 “ ~ ' (109) 

where ag = (ss) is the s quark condensate. The key ther¬ 
modynamic quantity, the partition function Z in three 
flavor case, is defined by 

Z{T, fij, fiB, gs, 1") = , (110) 


with baryon number, isospin number and strangeness 
number 


We examine now the possible pion and kaon superflu¬ 
idity at finite isospin and strangeness chemical potentials 
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(i^xs 7 °s , 


i?±(k) = Jk2 + M2±^, , 


(118) 



( 111 ) 


as conserved charges, and baryon, isospin and 
strangeness chemical potentials /xb,M 7 and fis- In the 
frame of imaginary time formulism of finite temperature 
field theory, the partition function can be represented as 


ZiT, fii, fiB, fJ‘S,V) = 

, ( 112 ) 


where /i = diag{fj,u, ga, Ms) is the chemical potential ma¬ 
trix in flavor space with the chemical potential for the 
strangeness quark 

Ms = ^ - MS • (113) 


A. Chiral Properties at Low /ij and gs 

We discuss first the three chiral condensates (j„, ad and 
as under the condition of low isospin and low strangeness 
chemical potentials, m/ < M/ ^md ms < Msi where M/ is 
the critical isospin chemical potential for pion condensa¬ 
tion and Ms Ihs critical strangeness chemical potential 
for kaon condensation. They will be determined in the 
following subsection. Performing the standard mean field 
approach and keeping only the linear terms in the meson 
fluctuations, we obtain the Lagrangian in the mean field 
approximation 

d^mf = + m 7 o - M)ip 

- 2G{al + al + a'^) - AKauadas , (114) 

where M = {Mu, Md, Mg) is the mass matrix in flavor 
space with the effective quark masses 


Mi = moi - 4,Gat + 2Kajak , 

{i = u,d,s, iydj^ k, j < k) . (115) 

Because of the lack of pion and kaon condensates, the 
mean field quark propagator is diagonal in flavor space. 


Smf{k) = diag(Su{k),Sd{k),Ss{k)'j 

with the matrix elements 
A+7o 


S^ik) = 


-f 


Ah 70 


ko-E-{k) ko + E+i'k)' 


(116) 


(117) 


where Ef^ are the effective quark energies 


and the energy projectors 

^ 7o(7-k-fM,) ^ 

^ Vk2 + Mf ) 


(119) 


In self-consistent Hatree-Fock approximation the gap 
equations which determine the value of the chiral con¬ 
densates au, ad, as are expressed in terms of the quark 
propagators, 

/ 74 

Jj^^'^DS^{p) . (120) 

After performing the Matsubara frequency summation, 
we have 


ai = -6 


d^k Mi 

{2n f y/p m2 


(1 -/(f;-) -/(£;+)) . 


( 121 ) 


The three effective quark masses are shown in Fig. (23) as 
functions of baryon chemical potential mb at T = /rs = 0 
and Pi = 0.06 GeV < m^r. We use the model parameters 
fixed in [46]. The three masses keep their vacuum values 
till the common critical point where they suddenly jump 
down. The mass difference between u and d quarks arisen 
from the finite isospin chemical potential can only be seen 
after the chiral transition. 



Hb (MeV) 

FIG. 23. The three effective quark masses M„, Md and Ms 
as functions of baryon chemical potential /is at T = /rs = 0 
and gi = 0.06 GeV in the flavor SU(3) NJL model. 


In principle the meson mass spectra can be evaluated 
by solving the corresponding pole equation like (54) for 
two flavor NJL model. The meson polarization functions 
in flavor space are evaluated in Appendix D. In gen¬ 
eral case with nonzero pion and kaon condensates, the 
off diagonal elements in the polarization function matrix 
in flavor space makes it complicated to solve the pole 
equation. However, in the region without pion and kaon 
condensates we are interested in here, the off diagonal 
elements disappear and the pole equation is reduced to 
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l-2G±nMM(fco =MM,k = 0) =0 , (122) 

with for M = a°, a'^, a~, for M = k~, Gq for 

K®, Gi for M = 7 r°, 7 r+, 7 r“, Gj for M = K~^, K~ and 
Gq for . At zero temperature, we can analytically 

obtain the meson masses as functions of /i/ and /xg by 
comparing the above mass equations at finite /x/ and /xg 
with the same equations but in the vacuum, 

fJ^s) = mM — QI^J-I — QsfJ'S , (123) 

where Qi and Qs are the isospin number and strangeness 
number of the meson M, respectively. 

B. Phase Diagram in /xs — /x/ plane 

In this subsection, we analytically determine the criti¬ 
cal isospin chemical potential /xj for pion condensation 
and the critical strangeness chemical potential /Xg for 
kaon condensation in the three flavor NJL model, and 
briefly discuss the phase structure in /xg — plane. For 
simplicity, we always set T = hb = 0 in this subsection. 

We first consider the case with only pion condensation 
in the region of |/xs| < /Xg and determine the critical 
value Taking the standard mean field approximation, 
we obtain the inverse of the quark propagator 

/ G-i 2 zG+75^ 0 \ 

S-\ik) = 2zG+757r G,-i 0 , (124) 

V 0 0 Gji J 

with G~^ = + jii'jo — Mi for i = u,d, s. Since in the 

current case the pion condensate is decoupled from the s 
quark propagator, we can take the same procedure as in 
the two flavor case and obtain the same critical isospin 
chemical potential 

/xf = (^s < /xg) . (125) 


At mean field level the inverse quark propagator in the 
flavor space can be expressed as 

/ G-i 2zG+i^„,75 0 \ 

= 2iGtK^sl^ Gji 0 . (128) 

V 0 0 Q' J 

For convenience, we have changed the basis in the flavor 
space from (zx, d, s) in (124) to (zx, s, d) in (128). It 
is straightforward to derive the propagator itself from its 
inverse, 

/ Su{k) Sus{k) 0 \ 

S^f{k) = Ssu{k) Ss{k) 0 , (129) 

V 0 0 Sd{k) ) 

with the elements 

Su{k) = + d«7o - Mu) (130) 

-f 4(G|)2A:2^75(7''fcM + MsTo - Ms)~^-1b 
S d{k) = + ^d7o - Md)~'^ , 

Ssik) = + ^s7o - Ms) 

+ 4{Gt)‘^Klsl5{l'"k^ + fiujo - Mu)~^J5 
Sus{k) = -2iG^KusSu{k)j5{'y^kf, + ^slo - Ms)~^ , 

Ssu{k) = -2iG^KusSs{k)"f^{'-f^k^ + /x„7o - M„)“^ . 
With the propagator we can express the condensates as 

a, = -Nc J ^^TTD[iiS^{k)] , (z = zx,d,s) 

Kus = Nc J [{Sus{k) + Ssuik)) 75 ] . (131) 

To calculate the critical strangeness chemical potential, 
we set Kus O'*', and the gap equations on the phase 
transition line become 


We then study the case with only kaon condensate in 
the region of < m^. We define the K~^ and K~ con¬ 
densates 

Kus = 2(ziz75s) = 2(sz75zx) , (126) 

and , K'^ condensates 

Kds = 2(Jz75s) = 2(sz75(i) . (127) 

Which kind of kaon condensation happens depends on 
the sign of /xg and /x/. In the physical world with /x/ < 0, 
only and K~ condensates can be realized. Without 
losing generality, we consider here the case /x/ > 0 and 
/xs > 0 only, and then map the result to the other regions 
in the — d/ plane. In this case the possible kaon 
condensation is 77+. 


(Ji -f iNc [ —rrTrc- = 0 , 

J (27r)4 - Mi 


(z = u,d,s) , 


l + 2iV,zG+ J 


--7=; 

+ M«7o - Mu 7+fc^ -I- Hsjo - Ms 


(132) 


= 0 . 


Doing the trace in Dirac space, performing the Matsub- 
ara frequency summation, and then comparing the gap 
equations with the mass equations for kaons in the vac¬ 
uum, we find 


fig = niK - fii/2 , [fii < m^) . (133) 


This result is certainly independent of the model param¬ 
eters and the regularization scheme. The above analytic 
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conclusions /ij = and fj,g = rriK — fJ^i were obtained 
numerically in the three flavor NJL model without Ua{^) 
breaking term [18]. 



-400 -200 0 200 400 

g, (MeV) 

FIG. 24. The phase structure of pion and kaon condensa¬ 
tions in jj,s — plane at T = hb = 0 in the flavor SU{3) 
NJL model. The upper panel is for fii > 0 and jj,s > 0 and 
the lower panel is for all possible jj,i and ps- 

Combining the two critical chemical potentials /i/ and 
/i 5 , the system is in the normal phase without pion and 
kaon condensations when fxj and /ig satisfy the con¬ 
straints 


<m^ , 

^J-S < rriK - IJiil‘2- . (134) 


In the region of /i/ < tUtt and /xg > rriK — /ij/S, 
there is only K'^ condensation, and in the region of 
/is < rriK — /‘^ and /i/ > niTr there is only pion conden¬ 

sation. The phase transition lines in the fis ~ IJ'i plane 
are indicated by solid lines in Fig.(24). As for the other 
region in the /is — plane we should consider both pion 
and kaon condensations and the calculation will become 
much complicated. The dashed line in Fig. (24) which 
separates the region with pion condensate from the re¬ 
gion with kaon condensate is just an estimation. How¬ 
ever, with the bosonized chiral Lagrangian, this dashed 
line is described by [12] 


Ms (Im/| > 





( 135 ) 


Considering the symmetries between 7r+ and 7r_ and 
among K'^, K~, and K^, we display in the lower 
panel of Fig. (24) the phase diagram with different pion 
and kaon condensations, which is the same as the one 
obtained in the frame of bosonized chiral lagrangian [12]. 


IX. SUMMARY 

It is well-known that physical symmetries of a system 
dominate its behavior not only in the vacuum but also in 
hot and dense medium. The symmetries in the vacuum 
are usually changed by the temperature and density ef¬ 
fect. In this paper we have investigated the changes of 
isospin symmetry, chiral symmetry and color symmetry 
and their reflection in meson spectra and thermodynamic 
functions at finite temperature T and isospin, baryon and 
strangeness chemical potentials /t/, /is and /xg in effective 
models. 

We studied the flavor SU{2) NJL model in detail. 
With the standard technics we set up in flavor space 
the quark propagator matrix with off diagonal elements 
arisen from the pion condensation at finite isospin den¬ 
sity. With the propagator we derived the gap equations 
determining the temperature and chemical potential de¬ 
pendence of the pion and chiral condensates, and ob¬ 
tained the thermodynamic functions in mean field ap¬ 
proximation. By self-consistently solving the gap equa¬ 
tions at finite isospin density and the pion mass equation 
in the vacuum, we analytically derived the critical isospin 
chemical potential for pion superfluidity, /xj = m^, at 
quark level. Therefore, in chiral limit with vanishing pion 
mass, any small isospin effect will lead to chiral symmetry 
restoration and isospin symmetry spontaneous breaking 
even in the vacuum. In real world with finite current 
quark mass and at finite temperature, there is a strong 
competition between the pion and chiral condensates in 
their coexistence region. The chiral condensate goes up 
continuously in the mixed region and drops down out¬ 
side, while the pion condensate decreases monotonously 
in the whole temperature region. The bulk properties of 
the system in the limit of high isospin density is also quite 
different from the corresponding ideal gas, the equation 
of state is hard and the ration of pressure to energy den¬ 
sity approaches 0.7. We investigated the effect of pion 
superfluidity on the meson properties in the hot and 
dense quark-meson plasma, by taking into account the 
off diagonal elements in the meson polarization function 
matrix in isospin space. We found that the mean field 
approximation to quarks together with RPA to mesons 
can well describe the isospin spontaneous breaking, there 
is a massless Goldstone mode which is a linear combina¬ 
tion of the 7r_|_, 7r_ and a modes in the whole superfluid¬ 
ity phase. We proved the Goldstone mode analytically in 
general case with finite T, fii and /is. In the superfluidity 
phase, while the ttq mode is still an eigen mode, the ct, 7r_|_ 
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and 7r_ mixing is very strong, especially in the region 
close to the phase transition point. We discussed also 
the competition between the pion superfluidity and color 
superconductivity at finite baryon density. While we did 
not make full calculation with both pion and diquark 
condensates, we determined the critical baryon chemical 
potential for color superconductivity. The region of 
the normal phase without pion and diquark condensates 
in —p./ plane is controlled by the pion mass ttItt and di¬ 

quark mass mu in the vacuum. In the case of two colors 
it is in a square box defined by /x/ < ttItt and /is < 
and in the case of three colors it is in a rectangle defined 
by Hi < and hb < l.Smc. 

We briefly extended our discussion on pion superflu¬ 
idity to the flavor SU{3) NJL model at finite isospin 
and strangeness densities. We found that the conclu¬ 
sion of the critical isospin chemical potential 
for pion superfluidity does not change, and that the crit¬ 
ical strangeness chemical potential for kaon superfluidity 
is related not only to the kaon mass rriK but also to 
the isospin chemical potential, Hs = ~ We 

showed the phase diagram of pion and kaon superfluidity 
in Hs — fJ'i plane with possible charged pion condensates 
and kaon condensates. 

We have investigated the pion superfluidity and me¬ 
son properties at tree level and the thermal excitation 
at finite isospin density in the bosonized version of the 
NJL model, namely, the linear sigma model, and made 
comparison of our calculations in the NJL model and 
sigma model with the result obtained in chiral perturba¬ 
tion theory. Many conclusions in the linear sigma model 
are qualitatively or even quantitatively consistent with 
that in the NJL model. We found that at zero tempera¬ 
ture the linear sigma model can reproduce the result of 
chiral perturbation theory by taking the limit ma- —>■ oo, 
but at finite temperature the Hartree-Fock approxima¬ 
tion can not recover the Goldstone mode in the pion su¬ 
perfluidity phase. The phase diagram in the T — hi plane 
in large N expansion approach is more reasonable than 
the one in the Hartree-Fock approximation. 

It is necessary to discuss the difference among the ef¬ 
fective models. The conclusion h‘i = is observed 
in almost all model calculations and lattice simulation. 
While obtaining this result is easy and looks trivial in the 
meson effective models such as chiral perturbation theory 
and linear sigma model, because the pion mass is a pa¬ 
rameter in these models, it is in the NJL model a general 
result of mean field approximation to quarks and RPA 
to mesons and diquarks. This is also true for kaon con¬ 
densation. In chiral perturbation theory, the chiral and 
pion condensates are regarded as a rotation from each 
other. As a result, the maximum pion condensate should 
be the same as the chiral condensate in the vacuum. In 
the NJL model, we found that this rotation point of view 
is only true around the phase transition point. Beyond 
this region the pion condensate can be larger than the 


maximal chiral condensate in the vacuum. This behav¬ 
ior is qualitatively consistent with the lattice simulation 
[9-11]. As for the linear sigma model, because the pion 
condensate increases with hi without bound, the phase 
structure at high hi is unphysical. Of course, due to 
the lack of confinement mechanism in the NJL model, 
we need a momentum cutoff A in numerical calculations. 
Therefore, we can not take the numerical results in high 
temperature and high density regions seriously. 
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APPENDIX A: TRACE IN DIRAC SPACE 


In obtaining the gap equations for the condensates and 
the meson and diquark polarization functions, we need to 
take the trace of the quark propagators including energy 
projectors in Dirac space. Here we fist the results. 


Tr[A+(p)7oA+(p)7o] 

Tr[A+(p)7oA_(p)7o] 

Tr[A+(p)75A+(p)75] 

Tr[A+(p)75A_(p)75] 

Tr[Ae(p)7oAe/(p)75] 

Tr [75A4 .(p)7o75A+(p)7o] 
Tr [75 A+(p)7o75A_(p)7o] 
Tr [75 A+(p)7575A+(p) 75] 
Tr [75 A+(p)7575A_(p) 75] 
Tr [75Ae(p)7o75Ae'(p)75] 
Tr[Ae(p)7o75Ae'(p)7o] 

Tr[Ae(p)7575Ae/(p)75] 

Tr[A+(p)7o75A+(p)75] 

Tr[A+(p)7o75A_(p)75] 

Tr[A_(p)7o75A+(p)75] 


= Tr[A_(p)7oA_(p)7o] 

= 2M‘^/El , 

= Tr[A_(p)7oA+(p)7o] 

= 2 - 2M'^/El , 

= Tr[A_(p)75A_(p)75] 

= 2 - 2M'^/El , 

= Tr[A_(p)75A+(p)75] 

= 2M‘^/El , 

= Tr[Ae(p)75Ae-(p)7o] 

= 0 , 

= Tr [75 A„(p)7o75A_(p)7o] 
= 0 , 

= Tr [75 A_(p)7o75A+(p)7o] 
= -2 , 

= Tr [75 A_(p)7575A_(p)75] 
= 2 , 

= Tr [75 A„(p)7575A+(p)75] 
= 0 , 

= Tr [75Ae(p)7575Ae'(p)7o] 

= 0 , 

= Tr [75Ae(p)7oAe/(p)7o] 

= 0 , 

= Tr [75 Ae(p) 75 Ae/ (p)75] 

= 0 , 

= Tr [A_(p)7o75A_(p)75] 

= 0 , 

= 2M/Ep , 

= -2M/Ep , 
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Tr [A+ (p) 7575A_ (p)7o] = Tr [A_ (p)7575A+ (p)7o] 

= 0 , 

Tr[A+(p) 7575 A+(p) 7 o] = 2 M/Ep , 
Tr[A_(p) 7575 A_(p) 7 o] = - 2 M/Ep , 

Tr [75 A+(p)75A+(p)7o] = Tr [75A_(p)75A_(p)7o] 

= 0 , 

Tr[ 75 A_(p) 75 A+(p) 7 o] = 2 M/Ep , 
Tr[ 75 A+(p) 75 A_(p) 7 o] = - 2 M/Ep , 

Tr [75 A_ (p)7o A+ (p)75] = Tr [75A+ (p)7oA_ (p)75] 

= 0 , 

Tr[ 75 A+(p) 7 oA+(p) 75 ] = 2 M/Ep , 
Tr[75A_(p)7oA_(p)75] =- 2 M/Ep. (Al) 


APPENDIX B: POLARIZATION FUNCTIONS IN 
SU(2) NJL MODEL 

We list the meson polarization functions 11^ (fc) in the 
SU(2) NJL model for i,j = ct, tti , 7 r 2 , tts and for i,j = 
CT, 7 r+, 7 r_, ttq. With respect to the basis (ct, tti, 712 , 713 ), 
they are 

noo(fc) = iNc J -^^TTD(Suu{q)Suu{p) + Sddiq)Sddip) 

T ^ud{q^^du{p} T ^du{q)^ud{p}^ ; 

Tliiik) = -iNc J -^^"^^D{lbSuu{q)l5Sdd{p) 

+ '1bSdd{q)lbSuu{p) + '1bSud{q)lbSud{p) 

+ lbSdu{q)i5Sdu{p)^ , 

Ti22{k) = -iNc J -^^TTD^jdSuuiqhbSddip) 

+ IbSddiqhdSuuip) - lbSud{q)lbSud{p) 

- lbSdu{q)l5Sdu{p)^ , 

n33(fc) = -*^c J -^;^TTD(j5Suuiqh5Suuip) 

+ 'YbSddiqhbSddip) - 'ybSudiqhsSduip) 

- 75Sdu{qh5Sud{p)j , 

Tloiik) = -Nc J ■^;^TTD(Suu{qh5Sduip) 

+ SddiqhbSudip) + Sud{q)ibSuu{p) 

+ Sdu{q)i5Sdd{p)^ , 

nio(A:) = -iVc J ■^;^T^T^D(j5Suu{q)Sud{p) 

+ 'r5Sddiq)Sduip) + i5Sud{q)Sddip) 
y5^du(q)^uu(p)'j ^ 


no2(k)=iJVc J -^^^TToiSuuiqhbSduip) 

- SddiqhbSudip) + SduiqhsSddip) 

n2o(fc)=*A^c J ■^i^T^^D{ibSdd{q)Sdu{p) 

- lbSuu{q)Sud{p) + lbSdu{q)Suu{p) 

- i5Sud{q)Sdd{p)^ , 

no3(fc) = iVc J -^^T^T:D{Sdd(.q)i5Sdd(p) 

- Suu{q)lbSuu{p) + Sud{q)lbSdu{p) 
— Sdu{ 


n3o(fc)=A^c J ^^^TrD(^755d<i(g)5dd(p) 

- 'y5Suuiq)Suuip) + i5Sduiq)Sudip) 
^d^ud{q)^duip)^ 1 

ni2(fc) = -iVc J ■^;^TTD(l5Suuiqh5Sdd{p) 

- 'y^SddiqhdSuuip) + ^bSduiqhsSduip) 

- 75Sud{q)l5Sudip)j , 

n2l(fc) = -iVc J -^;^TTD(l5Sdd{q)l5Suu{p) 

- 'ybSuuiqhsSddip) + 75Sduiqh5Sduip) 

- 75Sud{q)l5Sudip)j , 

ni3(fc)=iA^c J ^^^TrD(^755„u(9)755„d(p) 

- ibSdd^qhbSduip) + i5Sdu{q)i5Suu{p) 

- ibSud{q)ibSdd{p)) , 

n3i(fc)=*A^c J ^^^TrD(^755„u(9)755d„(p) 

- lbSdd{q)lbSud{p) + lbSud{q)lbSuu{p) 

- i5Sdu{q)i5Sdd{p)^ , 

/ d^ 

j^^"^^D{lbSdu{q)lbSuu{p) 

+ IbSudiqhbSddip) - lbSuu{q)lbSud{p) 
- JbSddiqhbSduip)^ , 

n32(fc)=iVc J ^^^TrD(^755„u(g)755d„(p) 

+ lbSdd{q)lbSud{p) - lbSud{q)lbSuu{p) 
- i5Sdu{q)i5Sdd{p)^ , 


(Bl) 


where the momentum q is defined as q = p k. 

For determining the meson masses in RPA, we need 
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only the polarization functions at k = 0. With the cal¬ 
culated trace shown in Appendix 1, we have 

no3 = Ilao = Ilia = Ilai = 1123 = 1132 = 0. (B2) 

If we further set ko = 0 which is used to determine 
the temperature and chemical potentials corresponding 
to the zero meson mass, we obtain 


U^^aik) =-V^Nc j Trg (jbSuu (p + k)Sud(j>) 

+ IbSudip + k)Sdd{p)j , 

U^_a{k) = -V2Nc J ■^^TTD(j5SduiP+ k)Suuip) 

+ IbSddip + k)Sdu{p)^ ■ (B5) 


Bn? = 1120 = ni2 = 1121 = 0. 


(B3) 


After performing the Matsubara frequency summation 
we derive the expressions as explicit functions of temper¬ 
ature and chemical potentials, 

n22(fco = 0) = -2iV, J {f{EZ) - 




(B4) 


The explicit expressions as functions of temperature and 
chemical potentials at k = 0 can be written as 




r{ko)=2N,J 


(fp 


(271)3 Ej 

E-E+ -{El-p.yA)-AG^n^ 
kl-{EZ-E+f 
(fiEZ) + f{-E+) - /(-A+) - fiEt)) 


E+ 


En 


Ilssiko) = 2N, J 


f dZp 

^ EZE+ + {El - pH A) + 4G2712 / 1 ^ 1 \ 

J (27i)3 

kl-{EZEE^f \E^ EZ J 


' E-E+ - {El - / J_ 

kl-{EZ-E+Y \Ei: EZ J 

(fiEZ) + fi-Et) - fi-EZ) - fiE+)) 

^ EzE+ + {El-t4/A)+AG^n^ / J_ ^ 

k^-{EZ+E+)^ \E^ EZ J 

(f{EZ) + f{E+)-f{-EZ)-f{-E+)) . 

The meson polarization functions in the basis 
(ct, 71+, 7r_, tto) are defined as 


U„„{k)-iNcj (-2^)4( ^UU {q)Suu{p) + Sddiq)Sddip) 
T ^ud{q^^du(^p) T ^du(^q)^ud(^p)^ ? 


n7j-_^7r4 

Xk) 

= -2iN^ J 

f d‘^P 

(271)4 

Tru (zj5Suu{p + k)^i,Sdd{p)^ 

n7r_ 7r_ 

Xk) 

= -2iNc J 

f dXp 

(271)4 

Tru {z5Sdd{p + k)"f5Suu{p)^ 

n7r^7r_ 

Xk) 

= -2iNc J 

f d'^P 

(271)4 

Tru (j5Sud{p + k)-f5Sud{p)^ 

n7r_ 774 

Xk) 

= -2iNc J 

f d'^p 

(271)4 

Tru (^75 

u{p + k)"f5Sdu{p)^ 

n(77r_|. 

Xk) 

= -V2Nc 

f d'^p 

J (27r) 

jTrz)(^5„^ 

{p + k)^5Sdu{p) 



+ Sdu {p + 

k)z5Sdd{p)^ , 


n^-TT- 

{k) 

= -V2Ne 

f d-p 

J (27r) 


{p + k)^5Suu{p) 


+ Sddip + k)"f5Sudip)^ , 


{f{EZ) + f{Et) - f{-EZ) - /( - E+)) 


+ANc 


d^p M2 


(271)3 El 


4G27i2 1 


kl - A{EZY Ep 


^(/(AI)-/(-£;;)) 


4G27i2 


kl - AiE+f E+ 


if (El) - f{-El)) 


n^+7i-+(A:o) - ANc J ^2^^3 

'(A-)^ + (Ep - + fco {Ep - III 12) 

kl-A{Ezf 

^if{EZ)-f{-EZ)) 

h/p 

{E+f + {Ep + ^ii/2f - fco {Ep + /ij/2) 
^ kl-A{E+f 


^(/(A^)-/(-£;+)) , 

H/p 

n„_(W = 4+||A 

'(A-)^ + {Ep - ^ii/2f - fco {Ep - III 12) 

kl-A{Ezf 

-^{!{EZ)-!{-EZ)) 

h/p 

{Ep^ + {Ep + Hi/2) A- ko {Ep A- pi 12) 

^ kl-A{E+Y 
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E+ 


n^+7i-_(A:o) = n^_7i-+(A:o) = -4iVc J 


d^p 


1 


lk^-4{Ep)^Ep 
4G27r2 1 


^(/(£;z)-/(-£+)) 

{f{E+) - f{-E+)) 


kl - A{E+Y E+ 

^a^+{ko) = n^+o-(fco) = 4v^iVcG7r J 

2Ep — Hj + fcg 1 


d^P Mq 
(27r)3 £p 


fc2-4(Sp-)2 
2£/p + /i/ — fcg 1 


^(/(£;z)-/(-i?+)) 
(/(i?i) - )) 


A:2-4(£;p+)2 E+ 

U^^Jko) = U^Mko) = ^V2 N,Gtt f 7 ^^ 

J [2tt) Ep 


and the effective potential 

m = - niDe+ 


d ' 
2 . 


a a) + 2 {'ktt) + 


+ 2( {aaY + 2 {'K'kY + (ttotto)^) 


(C3) 


(B6) 

where the thermal excitation functions {era), (tttt) and 
(ttotto) are calculated from the mean field propagators, 




1 


1 


{aa) = 


(27r)3 + M2 gVk^+M2/T _ ^ 

/ d^k 1 1 

(27r)3 y/p + ^ gVk^Z^/T _ J 
/■ d3k 1 
J (27r)3 2y/fc2 + m2 

1 1 


(C4) 


j(Vk^+Mf-A«j-)/T _ g(Vk=+M2+pj)/T _ ^ 


The thermodynamic potential of the system in Hartree- 
Fock approximation can be expressed as the summation 
of the quasi-particle contributions plus the effective po¬ 
tential, 


APPENDIX C: HATREE-FOCK 
APPROXIMATION IN LINEAR SIGMA MODEL 

In Hatree-Fock approximation, the interaction terms 
among a and tt^ in the effective Lagrangian (81) of the 
linear sigma model are fully absorbed into the a and tt^ 
effective masses, and the model behaviors like a quasi¬ 
particle system. At sufficiently high temperature, the 
pion condensation vanishes and the system contains only 
chiral condensate. This is the case we considered in de¬ 
riving the critical isospin chemical potential /xj for pion 
superfluidity in Section VI. In this case the Lagrangian 
density is simplified as 

Ehf = ^(^{da)'^ + ((dTTa)^ -h {dtni - gnT2f (Cl) 

+ (9(7X2 -f giTTif - [VlTif + (V7r2)^ ) 

- \ [M^a^ + M^ttI + M^iwl + 7 r|)] -U {0 ■ 


n = c7(^) -h T 



1 _ 

0; 




d^k 

(27x3) 

d3k 

(27x3) 

d3k 

W) 


In 


In 


In 


1 _ g-Vk"+Mg/T 


1 _ g-(VkG^-w)/T 


1 _ g-(Vk2+M2+(t,)/T 


(C5) 


The gap equation for the chiral condensate ^ is derived 
by minimizing the thermodynamic potential, 

25 ^? - /^ + 3(crcr) + 2(7X7 x) -h (txqtxo)) -f ^ml = f^ml . 

(C6) 


The mass equations (C2) with the thermal excitation 
functions (C4) and the gap equation (C6) form a group 
of coupled equations and determine self-consistently the 
temperature and chemical potential behavior of the chi¬ 
ral condensate and the effective meson masses. 


with the effective meson masses 

= ‘^gl - fl + 3(o-cr) + 2(7X7 x) -F (txqtxo)) -f ml , 
Mq = ‘2gl - fl + {aa) + 2(7X7x) -h 3(7Xo7Xo)) -f ml , 

= ‘^gl - fl + {^^) + 4(7!'7x) -h (txqtxo)) -f ml , 

(C2) 


APPENDIX D: POLARIZATION FUNCTIONS IN 
SU{3) NJL MODEL 

We define the following matrices 

A± =^(A,±zA,) . (Dl) 
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The meson polarization functions are defined as 

^MM{k)=i J -^^TT{rijSmfip + k)TMSm,f{p)) , 

(D2) 


with the couplings 


' A ^2 M = a+ 
Xi 2 M = a~ 
A 3 M = a° 
A+ M = K+ 
A 45 M = K~ 
A+ M = K° 

Ta^ = < "^67 M = R° 

*75^12 Af = 7r+ 
* 75 '^r 2 M = TT~ 
775 A 3 M = TT° 
*75A45 M = K+ 
*75A45 M = R- 
775A+ M = K^ 

^ 775A^7 M = K^ 

' X[2 A7 = a+ 
A+ M = a- 
A 3 M = 0 ° 
A 45 M = K+ 
A 45 M = K~ 
A -7 M = kO 
p* _ A^y M = 

^ 775 A ^2 A7 = 7r+ 

775 A+2 ^ = TT 
775 A 3 M = TT^ 
*75A45 M = K+ 
*75A45 M = K~ 
M = K° 

^ 775A+ M = ifO 


(D3) 


Doing the trace in color and flavor space first, we obtain 


^aoaoik)=iNc J -^^Ttd(Su{p +k)Su{p) 

+ Sdip + k)Sd{p)^ , 

T^a+a+ik) =2iNc J ■^^TTD(Suip + k)Sdip)j , 

na_o_(fc)=67 J -^^TTD(Sd{p + k)Suip)j , 

^k+k+{k) = 2Nci J ■^^TTD(Suip + k)Ssip)j , 

T^k.k.ik) =2Nci J -^^TTD(Ss{p + k)Suip)j , 
TlkOko{k) = 2Nci J ■^^TTD(Sdip + k)Ss{p)j , 
Il-,okoik)=2NJ J ^Trd,(Ss{p + k)Sd{p)) , 


^TTo^oik) = -Nci J ^^^TrD(^755„(p + fc)755„(p) 

+ IbSdip + k)j5Sd{p)^ , 

n^+7r+(fc) =-2iVci J ^;^Ttd{j5Su{p +k)-f5Sd{p)j , 

/ d^ 

(27r^4 Trji (j5Sd{p + k)^5Suip)^ , 

^K+K+{k) = -2Nci J ^^^TrD(^755„(p +fc)755s(p)) , 

/ d^ 

^2:^TrD(^755 s(p + fc)755„(p)^ , 

/ d^ 

(Ts^db + fc)755s(p)) , 

^^;^TrD (-f5Ss{p + k)-fr,Sd{p)j ■ 

(D4) 
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